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ABSTRACT 


A modern computer based automatic data acquisition/con- 
trol system was installed at the Department of Aeronautics' 
Combustion Laboratory. This system utilizes an HP-85 desktop 
computer as system controller for the HP-3054A data acquisi- 
tion system. These major control applications include: 

(1) particle size determination by measurement of scattered 
laser light, (2) a vitiated air heater system, and (3) a 
solid fuel ramjet test unit. Details involving specific data 
acquisition/control techniques for these applications are 


seeseussed. 
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This thesis project was undertaken to develop and install 
a modern, computer based experiment control/data acquisition 
system at the Naval Postgraduate School Aeronautics Devart- 
ment, Combustion Laboratory. Previously, all laboratory 
experiments were manually controlled and all test data were 
reduced manually by hand with no real time data processing 
possible. Precise timing of test equipment and control 
features were difficult or impossible. 

The introduction of the HP-3054A Automatic Data Acquisi- 
tion system and the integration of all associated test equip- 
ment to a desktop computer have revolutionized current 
operations at the Combustion Laboratory. 

Large quantities of high quality data can be acquired at 
high data rates and stored for future use or processed real 
time in closed loop feedback control circuits for precise 
process control and display. 

Other electromechanical equipment such as Visicorders 
provide analog real time data and run in parallel with digi- 
tal equioment but can be remotely controlled for vorecise 
timing of experiments. 

Interactive programs on the HP-85 computer allow the ex- 
perimenter to do preliminary calculations, initial experiment 


set-up, and to perform accurate process control of system 
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functions. On line data reduction and display of results is 
performed by the system CRT display, the 7225B Plotter, or 
the system printer. 

Four basic experiments/experimental apparatuses are cur- 
rently automated using the automatic data acquisition system; 
a light scattering experiment, a vitiated air heater, a solid 
fuel ramjet, and a turbojet combustor test rig. The first 
three applications are discussed in subsequent sections along 


with the measurement techniques used to implement the process. 
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It. SYSTEM AND COMPONENT DESCRIPTIONS 


To produce accurate test results involving high speed 
machinery and processes and simultaneously provide precise 
timing and control signals, a modern data acquisition system 
has been installed at the Combustion Laboratory at the Naval 
Postgraduate School. 

Figure 1 shows a system block diagram of the Hewlett- 
Packard 3054A Digital Data Acquisition System. The heart of 
this system is the HP-85 microcomputer functioning as system 
controller, although any computer compatible with the Hewlett- 
Packard Interface Bus (HPIB) can be substituted or connected 
to the system bus. 

A description of each equipment connected to the bus is 


described below. 


A. 3054A DATA ACQUISITION AND CONTROL SYSTEM 

The 3054A is a computer-based automatic data acquisition 
and control system. The system 1s complete with a system 
mainframe, the HP-3497A, and two voltmeters (the HP-3456A 
and HP-3437A), that are interfaced via the 82937A HP-IB I/O 
card to the powerful HP desk top computer, the 85A. 

The 3054A System documentation has been designed to co- 
ordinate efficient use of each of the system instruments, 


and it includes a specialized software support package that 
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Figure 1. Data Acquisition System 





is designed to use the System instruments to achieve specific 


functions and operations. 


B. HP-3456A DIGITAL VOLTMETER 

The 3456A Digital Voltmeter is a 3% to 6% digit, inte- 
grating voltmeter with relatively high speed, 10 parts per 
million (PPM) basic accuracy and 100 nanovolt sensitivity. 
Its DC and True RMS AC voltage maximum of 1000 volts and re- 
Sistance measurement capability reduce the amount of signal 
conditioning necessary. It can detect 100 nanovolt changes 
in 100 millivolt signals at speeds of 48 readings per second. 
This capability is required for measuring thermocouples with 
the 3054A system to better than 0.1°C resolution. By se- 
lecting 4% digits of resolution, reading rate is increased to 
330 readings per second. The 3456A has selectable integra- 
tion times from 0.01 to 100 power line cycles and provides 
input guarding to give 140 db of common mode rejection, pro- 
viding for accuracy in the presence of noise. Two other 
noise reduction features are provided by digital averaging 
and an analog input filter. 

Up to 350 readings can be stored in the 3456A Read/Write 
memory in addition to providing a programmed time delay be- 
tween readings. A "voltmeter complete” signal provides the 
necessary means to synchronize the operation of the 3456A 
with the 3497A to allow analog scanning rates up to 300 


channels per second. 
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ee 35437A DIGITAL VOLTMETER 

The 3437A is a high speed, 3% digit, DC voltmeter which 
provides precisely timed sample and hold readings. It can 
perform a variety of tasks from scanning DC inputs to AC 
waveform analysis. Repetitive signals with frequency com- 
ponents up to 1 MHz and low frequency transients longer than 
500 usec can be rapidly digitized and analyzed. Parameters 
such as RMS value, residual DC, harmonic content, and peak 
woes can be readily and automatically determined. 

Using the HP-85 controller and a "Fast Handshake" trans- 
fer mode of operation, reading rates of up to 4000 readings 


per second are possible in the packed binary mode of operation. 


Bee o¢9/7A DATA ACQUISITION/CONTROL UNIT 
The 3497 data acquisition/control unit combines precision 
measurement capability with control functions. It can provide 
precision measurements of strain gauge outputs, thermocouples, 
pressure transducers, and other sensors, and can also provide 
Siegeea! interfacing and control. A digital clock is built 
into the 3497A mainframe for real time control applications. 
Five slots are available for plug-in I/O card options on 
the 3497A mainframe. Additionally, the 3497A capabilities 
may be expanded through the addition of one or more 3498A Ex- 
tenders. Up to 1000 analog channels and 1300 digital channels 
Can be obtained with the addition of 14 3498A Extenders. Each 


extender has 10 slots for plug-in I/O card options. 


te 





The 3497A mainframe also provides several key features 
that make it very suitable for use in computer controlled 
data acquisition systems. The following TTL compatible input 
and output terminals are standard on the 3497A: external 
trigger input, external increment input, break before make 
(BBM) sync pulse (input/output), voltmeter complete output 
signal, channel closed output signal, and a timer signal. 
Additionally, the 3497A has multiple programming commands 


providing versatility in process control applications. 


Be 3498A EXTENDER UNIT 

The capacity of 3497A mainframe may be expanded through 
the addition of from one to fourteen 3498A extender units, 
each having 10 slots in its frame for I/O modules of various 
types to meet specific application requirements. 

Special attention to the 3498A extender operating manual 
is required when adding additional I/O cards. Slot numbers 


peenrough 9 do not exist for digital cards. 


feeeene=o2901M FLEXIBLE DISC DRIVE 

The mass storage requirements for the system are handled 
by the 82901M dual-drive disc system which supports two 5% 
inch flexible double-sided, double density discs, providing 
a total of approximately 550K bytes. Drive #0 is established 
as the default mass storage medium. Over 220K bytes of mass 
storage are also available in the HP-85 internal casette tape 


@rive. 
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eee He—/7225B PLOTTER 

The HP-7225B plotter contains a 17601A personality module 
and is therefore a microprocessor based HP-IB plotter that 
mproe@uces high quality graphic plots on any size chart up to 
meeorx 297 mm. 

Pen movement is as small as .032 mm for high resolution 
plotting. The plotter instruction set includes 39 different 
Miser actions to provide such capabilities as point digitizing, 
labeling, character sizing, scaling, and window plotting. The 
7225 is interfaced through the HP-IB providing extreme 
Metrsatility. 

Pen velocity is programmable from 10 mm/sec to 250 mm/sec. 
Seven different dashed-line formats, symbol mode plotting and 


meer aqefined characters aid in trace identification. 


H. HP=-2631B LINE PRINTER 

The 2631B is a 180 character per second printer that pro- 
duces high resolution dot matrix printing. Characters are 
Printed from one of two 128 ASCII character sets. A 16 chan- 
nel fixed or programmable vertical forms control (VFC) is 
available for any desired page and text length with several 
print pitches to select from. An underline mode is also 
available. 

The 2631B uses a bidirectional printing head and, addi- 
tionally, monitors incoming data to determine optimum print 


Seeeetion for the next line. Other features include blank 


Le 





Space detection and suppression, eight print densities with 
underlining, automatic vertical and horizontal tab settings 
and margins. It also includes an end of line (EOL) wrap- 

around feature. When the printer receives data which would 
cause printing to exceed the right margin, it automatically 
moves the data (which would ordinarily be lost) to the next 


line. 


i. HP~85A COMPUTER 

The heart of the Data Acquisition system is the 85A desk- 
top computer which functions as system controller. The com- 
puter is interfaced with system components via the 82937 
HP-IB I/O card. The 85A uses enhanced BASIC programming 
language and supports the PLOTTER/PRINTER ROM, INPUT/OUTPUT 
ROM, MATRIX ROM, and the MASS STORAGE ROM, providing an ex- 
tremely versatile instruction set. The 85A also features a 
Seenlo CRT display with full graphics capability, an internal 
32 column thermal printer, live keyboard, buffered I/O, pri- 
Ority intercept, high speed tape cartridge, autostart, and 
eo trapping. 

It is also a versitile calculator and can be used with 
virtually no previous HP-85 experience. The user can halt an 
executing program using the pause key, perform a series of 
calculations on the keyboard, print out the answer, and then 
resume program execution by pressing the CONTINUE key. 

Reference 1 describes the HP-85 system characteristics and 


the extensive instruction set. 


ES 





J. HP-6942A MULTIPROGRAMMER 

For high speed data acquisition applications the HP-3054A 
is augmented with the 6942A multiprogrammer. The 6942A is 
capable of functioning independently from the 3054A as a 
Gaieter control unit for bidirectional communication, data 
processing, and data transfer between an HPIB controller (the 
HP-85 computer) and a multiprogrammer system. It can be used 
in a single-unit system employing from one to sixteen plug-in 
I/O cards, or in a multi-unit system consisting of one 6942A 
master unit and up to seven 6943A extender units. Each ex- 
tender unit can accommodate up to sixteen input/output (I/0) 
cards, allowing a multiprogrammer system to be expanded to 
mee 1/0 cards. 

The 6942A emplovs 32 different instructions in its in- 
struction set, permitting a wide variety of I/O card func- 
Etens Under program control of an HPIB controller. The 
Miiceamuiction set includes system control, output, input, card 
Semmeol, and system timing functions. System control in- 
structions establish the basic operating modes of the multi- 
programmer system. Output instructions are used to send data 
to output cards while input instructions are used to obtain 
Maeda ErLOM input cards. The card control instructions allow 
access to various circuits on the I/O cards to enable close 
control and monitoring of their various operations. The sys- 
tem timing instructions permit the sequencing of events and 


the measuring of elapsed time. A detailed description of the 
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instruction set is described in the 6942A User's Guide [Ref. 
2] as well as programming details for the HP-85 controller. 

In addition to the 16 I/O cards, an HP-IB interface board 
and a transmission system board are installed in the rear of 
the 6942A unit. The 6942A's internal mainframe circuitry in- 
cludes a CPU/ROM board and a control/RAM board. All power 
supplies for the I/O card slots and mainframe circuitry are 
built into the 6942A unit. The CPU/ROM board includes a 16 
bit microprocessor chip which decodes and executes all in- 
fee cions, and controls all data transfers. The 12K x16 bit 
ROM (Read Only Memory) contains the microprocessor controlled 
programs (firmware) required to process all instructions. 

The ROM also contains diagnostic programs for built-in self 
test and signature analysis/troubleshooting capabilities. 
The control/RAM board contains a 2K x16 bit RAM (Random Ac- 
cess Memory) which is used for the temporary storage of in- 
structions and data. 

Backplane control circuits on the control/RAM board de- 
code I/O card addresses and control I/O card functions. The 
HP-IB interface board provides the bidirectional communica- 
tions interface between the 6942A and an HP-IB controller. 

Three I/O card types are presently installed in the 6942A 
multiprogrammer: a 69736A Timer/Pacer card, two 69751A high 
speed analog to digital converter cards, and two 69790B 4Kx 


16 bit Random Access Memory cards. 
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1. 69736A Timer/Pacer Card 

At power-on the timer/pacer "wakes up" in the one 
shot mode which produces a single output pulse with a 
ecrystal-controlled duration programmable from one micro- 
second to about eighteen hours (65535 seconds). Another 
mode (the continuous mode) can also be programmed in which 
the output goes high for the programmed time interval, then 
low for the same interval, and continuously repeats this 
cycle until the card is reprogrammed or disabled externally. 
The card's single output pulse can be used aS a programmable 
delay or as an enable signal for a frequency counter card 
(69775A) in frequency measurement applications. The square 
wave output generated in the continuous mode can be used to 
pace A/D or D/A converters as they analyze or generate wave- 
forms or to drive frequency-to-voltage converters or other 
devices which require a programmable-frequency square wave 
mnput . 

Figure 2 [Ref. 2] shows a detailed block diagram of 
the timer pacer card. In the Light Scattering experiment 
the timer pacer card 1S programmed to the continuous mode 
and drives the external trigger input of 2 high speed A/D 
converter cards. The timer pacer is itself externally trig- 
gered by the Reticon Photodiode array #1 which also syn- 


chronizes the operation of Photodiode array #2. 
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wee ANalog tO Digital Gonverter Card #69751A 
Mae 69751A is a 12 bit, analog to digital (A/D) con- 


verter used to measure bipolar DC voltages in one of four 
momoges: +100 MV, t1 V, £10 V, or £100 V. A voltage measure- 
ment can be initiated by a programmed instruction or by an 
external trigger signal. Three manual range switches on the 
card select either the +100 MV, +1 V, or the +10 V range. 
Input voltages in the +100 V range are connected to the card's 
divide-by-ten attenuator. Sixteen optional DC. input ranges 

(8 bipolar and 8 unipolar) are also available. Instructions 
for modifying and recalibrating the A/D card for the optional 
DC input ranges are described in section 11 of Ref. 2. 

The A/D card uses a high performance sample and hold 
amplifier and a successive approximation A/D converter which 
provide a high conversion speed and excellent accuracy charac- 
teristics. The DC input voltage is guarded to maintain high 
input impedance. Optically coupled isolators are used to 
isolate the DC input voltage from the data lines. 

The 12 bit digital word, indicating the magnitude and 
Sign of the measured voltage, is stored on the card and is 
always available for readback to the multiprogrammer main- 
frame memory using card subaddress 9. As many as 1339 read- 
ings can be stored at one time in mainframe memory. The 
multiprogrammer firmware converts the digital word decimal 


form for readback to the controller. The controller can 
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read back the voltage value using appropriate HP-IB extended 
talk address. 

The use of Extended Talk Addresses and other second- 
ary addressing with the HP-85 computer requires the custom 
HP-IB command "SEND" which is described in detail in the 
HP-IB I/O programming manual [Ref. 3]. 

The 12 bit data word (voltage reading) is routed to 
the A/D's edge connector so that the A/D can be used in con- 
junction with the 69790B memory card to make measurements at 
rates that exceed the throughput of the controller. Up to 
33000 readings per second can be transferred from the A/D 
converter to the memory card via the edge connectors on each 
card. Each reading can be initiated by program commands or 
Mean External Trigger (EXT) signal. 

Either the 69736A Timer Pacer or the 69735A Pulse 
fiero Output card can be used to externally trigger the A/D 
converter. The 69736A Timer/Pacer 1s programmed in a con- 
venient seconds-milliseconds-microseconds format for pro- 
grammable delay or continuous triggering. 

Figure 3 shows the functional block diagram of the 
69751A Analog to Digital converter card. Various subad- 
dresses can also be written using the software described in 
section 5 of the User's Guide [Ref. 2]. 

3. Memory Card 469790B 
The 69790B Memory Card is a FIFO (first in first out) 


Perot that can handle a continuous stream of data that is 


24 








[7 




















Pats) ce icy Nb lelp) TeATHTq OF HoTeuy peeds yhtH--AowuerborAdtatnNwW “¢E sanhty 
NORIO] — 
ae CIA d 
ie 193138 ALYY Wd YI0NYL : 
a ASN@ 3 
SS3004d JO ON 71 
2033 7 
NOISH3ANOD 40 ON9 i 
es 
baa ) 
HPOOINL TWNHILX9 
f\ 
mwOO 1 NOWMINOD INGM SOTWNY O3L¥ 10S! 
O4WND IH s 
OOTWNY c= ua4ang |assanovans \| 2 
Non IM vie /y | 
ONIN In#c ; 
INO NILIVO 2 
N 
34408 | Aassauaqwans | y 
| VIHA Hy 
| 
| 
NENILIYO 
singino {2 
¥ivd g 
watts 
= tM SURAT O O/W YOI 
aalW ios =a Y3MOd G3L¥ 10S! 
TNH 2m <a 





se 1@” Sandans 
a31¥ 2051 


TYAS LMA 


SMe K-AtOO8eA SSwe VNaeoxwaid Zw 


iy 


a 





being transferred from the outside world to the controller 
or from the controller to the outside world. In another 
operating mode, the recirculating input mode, it can act as 
memeireular buffer that continuously stores data coming in 
from the outside world until it is programmed to stop ac- 
cepting new data. At this time it contains the last ‘n' 
words of data received, 'n' being the capacity of the card. 
One more operating mode, the recirculating output mode, al- 
lows all or part of the memory contents previously acquired 
to be read repeatedly by an external device. 

Each 69790B memory card consists of two plug-in 
cards that occupy two adjacent slots in the mainframe. The 
card designated card #1 occupies the left-hand slot and com- 
municates with card #2 through a cable. 

A memory card assembly can communicate bidirection- 
ally with the controller and with an external device con- 
nected to the edge connector on card #1 of the pair. External 
data transfers through 16 input or 16 output lines and are 
controlled by means of three input or three output handshake 
lines. Special purpose Memory Input (MI) and Memory Output 
(MO) instructions allow data transfers to be made between 
the controller and a memory card assembly without storing 
the data in mainframe memory. 

In a system that contains more than one memory card 
assembly, each assembly can handle data transfers to or from 


the outside world simultaneously with external data transfers 
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being made by the other assemblies. Absolutely no inter- 
actions occur between their external interfaces. An external 
interface is independent not only of the other external in- 
terfaces, but also of any other card or mainframe activity. 
Thus, within a single assembly, an external device can write 
to a memory card while the controller is reading from another 
part of the same card, or an external device can read from a 
memory card while the controller is writing to another part 
Sietme Same card. Through the use of controller interrupts 
or another means of program control, data can be loaded into 
memory and read from memory on a continuous basis in either 
Saeiiput Oran output operation. 

Figure 4 is a functional block diagram of the 69790B 


memory card. 
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Figure 4. Multiprogrammer--Memory Board [Ref. 2] 
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mers FARTIC&i® SIZE DETERMINATION BY LIGHT SCATTERING 
OF A HELIOM NEON LASER BEAM 


The purpose of this experiment is to determine the volume 
to surface mean diameters (D55) Or earerculate matter, in “the 
exhaust plume and combustion chamber of a solid propellant 
meeket motor. 

The equipment test set up is shown in Figure 5. Refer- 
ence 4 provides extensive information on the physical princi- 
ples involved and a description of the optics system. 

The helium neon laser is remotely controlled under soft- 
ware control and radiates a beam at a wavelength of 0.6328 
micrometers. The beam is then expanded and columated. It 
is then passed through a limiting orifice to reduce the beam 
diameter to 2.2 millimeters. 

The beam is subsequently passed into a beam splitter 
where one beam is turned through 90° and routed through a 
window in the rocket motor casing. The other beam is aligned 
directly in the exhaust stream of the rocket motor approxi- 
mately two inches behind the nozzle throat. 

Each beam is passed through a focusing lens and directed 
through a narrow pass filter to the edge of a linear array 
of photodiodes. 

Two EG&G RETICON G Series Solid State Line Scanners are 


used, each containing 1024 photodiode elements. The Reticon 


Zo 








weibetq waqyskgs uoTtzeutTuasjag SZTS eToTIAReRgd 


JOYLNOD JYvMOuWH 


TOYLNOD 
JYVMLAOS SidH— — 





"G SAnb rg 








| SSvd MO7 oon 


PHYAL MS 
SVd MO] YOLOW 


\ 
INV113d0Ud amos \\ 

ino o3aa AUS N 
}-—£ 








¢ +t AVYUV 
JIQOIGOLOHd 









NDING 1d 


YOUUIWN 


Y3aSvT 
®N -9H 


Y3ILLINdS WV3eE 


30 





RL1024G series photodiode array 1S optimized for solid state 
image sensor applications. It contains 1024 diodes arranged 
in a linear array on 25 micron centers. Diode #1 is normally 
set 1.54 millimeters off of the laser beam centerline to re- 
duce the amount of transmitted light present, and to allow 
for increased sensitivity to the diffracted light pattern. 

The G series array is connected through an RC-100/106 
Bemple and hold "Boxcar" circuit which provides its own 
Meek Oscillator for self scanning, although provisions for 
externally pacing the array with an external clock signal are 
installed. 

Photodiode array #1 (Fig. 5) is operated in the internal 
clock/internal start mode and is completely independent. The 
mmeetienal Clock oscillator 1s set to approximately 30 KHz. To 
prevent integrated "dark current" from making a significant 
contribution to the output charge the total time between 
start pulses should be kept less than about 40 milliseconds. 

mime) COUNtdoOwn circuit on the array board counts ‘'n' 
clock pulses to generate a start pulse in the internal start 


mode. The countdown switches are initially set according to 


the formula: 


L E 
where ty = time between start pulses in seconds 
f. = clock frequency in Hz 


switch setting 


Sal 





3 


Bio (or >): 30 (107) 


eng. n 


105.0 


n 
A countdown switch setting of 1080 to 1 will provide a time 
between start pulses of 36 milliseconds at a clock frequency 
of 30 KHz, and excessive “dark current" error build-up is 
avoided. 

Photodiode array #2 monitors the second laser beam which 
emerges from a window in the rocket motor body. Array #2 
Operates in the external clock/external start mode and is a 
Slave to photodiode array #1. Therefore, the countdown cir- 
cuits in array #2 are not used. Array #2 has only one out- 
put, the video "boxcar" signal representing the photodiode 
voltage output which is routed directly to A/D converter #2 
in the 6942A multiprogrammer. 

The boxcar voltage outputs from the photodiode arrays 
are routed through analog low pass filters which pass all 
Signal components less than 3000 Hz. 20 db gain is attained 
in the low pass filter amplifier. The filtered outputs are 
Boucea to the voltage input terminals (+W, -Y) of the high 
Speed analog to digital converters located in the 6942A 
multiprogrammer. When a series of readings are desired the 
HP-85 programs the 69736A timer to output a square wave output 
at approximately 30 KHz when externally triggered by the 
blanking pulse output from photodiode array #1. 

The square wave output from the timer pacer card is di- 


rected to the external trigger inputs of both 69751A high 
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speed (A/D) cards which begin making 1024 conversions of the 
analog voltage present at terminals (+W, and -Y). The 12 bit 
digital outputs from terminals A through N are hard wired to 
he Corresponding edge connector pins of the 69790B, 4K memo- 
ry cards associated with each (A/D) card. The memory is pro- 
grammed to accept a 12 bit format and to store 1024 readings 
fea first in first out (FIFO) input mode. The memory card 
is then armed to interrupt the computer after 1024 readings 
have been stored on the card. The (A/D) card performs a con- 
version of the analog voltage at its input terminals (W, Y) 
every time it observes a trigger from the timer pacer card at 
its external trigger terminals (m) at a rate of 30 KHz. The 
(A/D) card signals the completion of each reading by output- 
ting an end of conversion (EOC) and its complement (EOC). 

The EOC is connected to the data available (DAV) terminal of 
each 69790B memory board. 

The procedure employed for an experiment is outlined be- 
low. Figure 6 presents a flow diagram of the system. After 
applying power, the HP-85 is loaded with a program disc in 
drive #0 and a data disc in drive #1. The program named 
"HELP" is loaded and run to initialize variables and provide 
Meemmentation. At the completion of "HELP" the next program 
1s automatically loaded and executed by a process called 
chaining, which allows a program of virtually unlimited length 


to be run in the limited memory space of 32K of the HP-85 
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Mere 6. Light Scattering Program Flow Chart--Part One 





computer. "ACQDTA" now performs the process control, system 
timing, and data acquisition tasks. 

First, the laser is energized, data files are created 
with user specified file names, and initial readings are 
taken from both photodiode arrays by software command. Each 
array Simultaneously sends its output to its 69751A (A/D) 
converter in the 6942A multiprogrammer. Each (A/D) card is 
externally triggered every 33 microseconds by the 69736A 
Timer/Pacer card and converts the analog voltage present 
at terminals (+W, -Y). The resulting 12 bit digital out- 
Sijemis sent to a 69790B (4K) memory card. Each memory card 
is programmed to store sequentially 1024 readings and to 
interrupt the computer after 1024 readings have been stored. 
At this time the numerical data are transferred via the 
HP-IB to the controller, first array #1 then array #2. This 
raw data represents the "no particle" or "zero scattering" 
eenaition. 

The "no particle" numerical arrays are then stored as 
a data file on the disc to preclude loss of the data. 

ime computer then performs a continuity check of the 
firing circuit using the 3456A in a resistance measurement 
mode. After a successful continuity check, the 3456A is pro- 
grammed to record DC voltage. It is connected to the firing 
Switch while displaying (on the HP-85), "STANDING BY FOR 
IGNITION", and energizes the nitrogen purge for the motor 


windows. When the operator closes the ignition switch a 
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voltage greater than 10 volts is sensed by the 3456A and 
causes the computer to energize the drive motor of the visi- 
corder, which records the chamber pressure-time trace. 

Once the computer senses the ignition sequence it con- 
nects the 3456A to the pressure transducer in parallel with 
the visicorder trace. 

When the combustion pressure rises to a user specified 
threshold pressure in (PSI) the computer delays or waits a 
user specified number of seconds, and triggers the multipro- 
grammer to acquire a new data set from the photodiode arrays. 
When the data are taken a pulse output voltage (TTL) is sent 
to another visicorder channel as an event marker signal to 
indicate data acquisition time on the pressure-time trace. 
Readings taken during the rocket motor burn time will be 
different than pretest readings Since particulate matter is 
present in the gas stream. This data is transferred to the 
computer and then to a data file on the data disc as before. 

After 5 seconds the visicorder 1s disengaged, the window 
nitrogen purge turned off, and the laser is deenergized. The 
next program is chained. 

"CMBINE" is automatically loaded and run. This program 
executes disc operations which read the "no particle" data 
and the "particle" data into the computer. It then separates 
this data into photodiode array #1 data and photodiode array 
#2 data with the same file names as previously specified by 


the user. Each file contains 1024 records, each containing 
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two numbers, a "no particle" and a "particle" voltage. The 
computer prints a reminder of which data 1s in which data 
file and then chains the next program segment. 

“PLTDTA" plots the data from a data file specified by 
the user, plotting diode voltage output versus diode number 
(Fig. 7). The operator can elect to plot the other data 
mee Or load and run the data reduction program. 

The operator 1s then prompted to load and run "REDUCE" to 
begin the reduction of the raw data. Figure 8 shows a flow 
chart for the data reduction process. "Reduce" functions to 
compute a normalized intensity vector and an angle theta 
vector. The theta vector 1s computed based on diode number 
(1 to 1024), lens focal length, and the index of refraction 


of the media containing the particles by the formula: 


Ky + n(.025) 
§ = samara radians 
where K. = the distance from the centerline of the laser beam 


: to the centerline of diode #1 (mm) 


nm = diode number (1 to 1024) 
F = lens focal length (mm) 
me im@ex of refraction 
antes = 0 
water = 1.35 
.025 = distance between diode centers (mm) 
Upon examination of the "PLTDTA" results, the location 


of the peak intensities are printed. The user is then prompt- 


ed for the diode locations where the slope is to be determined 
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and then projected back to obtain a centerline value for the 
scattered laser power using a linear regression technique. 
First, the "zero scattering” data 1s subtracted from the 
"with particle" data. The resulting 1024 element vector is 
then regressed against diode number to obtain the intensity- 
axis intercept representing the centerline value by the fol- 


lowing formula: 


nA + BX 


it 
KG 


1 


ll 
KG 


XA + BX. D 


where A = intensity-axis intercept 


n = number of data points 
B = slope of intensity vs theta line 
n 
= ) 6. 
i i=1 1 
n 
z 2 
<a 
1=1 
} 
Y, = i: 
i i=1 1 
n 
Y, = L I.6, 
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and subsequently, a normalized intensity vector 1s formed as: 


: mer) = Yenp'*)) 
°NORM ' eG 
(Youp CLyp ) 
where i = 1 to 1024 
1 8NORM = normalized intensities 
Vey (1) = voltage with particles (vector) 
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Venyp (2) = voltage without particles (vector) 
eLp - Wennp) — o (previous calculation) 

The normalized intensity is then plotted on the 7225B 
plotter versus the theta vector in radians (Fig. 9). This 
data is then stored on disc for future use in a user speci- 
fied file. The operator is then asked if this run is to be 
a calibration curve; if yes, then a curve fitting program 
"CALIB" is automatically loaded and executed. 

"CALIB" performs a tenth order polynomial eater faGee EO 
the data set specified by the user in order to generate a 
smooth curve for future comparison purposes with both the 
universal (1.e. theoretical [Ref.4]) curve ("UNICUR") and 
other raw data runs. The user first enters the "D bar" 
value, i.e. the average diameter of the particles being used 
for the calibration. 


A theta bar vector is then computed from the theta vector 


as follows: 


» ie 
ay a a2 el 
a 
where ay = theta vector (1=-1024) (radians) 
D3 = D bar (microns) mean volume to surface diameter 


» 


.6328 (microns) wavelength 
Normalized intensity is then plotted versus the non- 
mimensional theta bar (9) as in figure 10. The user is then 


asked to enter the minimum and the maximum theta bar of 
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interest. The data is regressed to form a smooth curve fit 
with a polynomial of up to tenth order as specified by the 


user and taking the form: 


_ 2 3 
oe a ceonce | Doe +. . . + bee 


where the coefficients are obtained by solving the system of 


m+l1 equations 


} pa 1 am= | 
mo, + b CS Q 7 + b oi ne 
0 1j2) 1 252) a m;2, 1 ey 2 
n n n E n n 
2 m+ 1 

pee. + b, ) 6.° +b, } 8. + +b J e@. = j 1.9 

a1 + Bist > ean * = 5 

m mt+1 Zi > m 
oe. 6+ «#8 + : ee on 


Once the polynomial coefficients have been determined, a 
smootn 100 element curve 1s computed, plotted and stored on 
Maerailsc in a user specified file. 

The smooth curve is plotted on top of the intensity vs 
theta bar plot (Fig. 10) to show the correlation, which is 
routinely excellent. If desired, this data file containing 
smooth curve data can be saved permanently on the floppy disc. 

"REDUC2" is automatically loaded and executed upon com- 
mereron Of "CALIB". The function of “REDUC2"” is to plot the 


meer sal theoretical curve (datafile "UNICUR") on a three 
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Sycle semi-logarithmic scale (Fig. 11). After the theoretical 
Suave 1S plotted the user may elect to plot any other calibra- 
tion curve for comparison purposes (Fig. 11). 

At the conclusion of "REDUCE", if the operator had not 
designated a run as a calibration run, the next data reduc- 
tion program is chained. 

meoouUCS” functions to compute D235 from the normalized 
intensity and theta data by comparing the data to specified 
reference curves (UNICUR, etc.) at 15 preset intensity levels. 
An average D.. is computed across this range of intensity 
levels. REDUC3 has tabular output only. 

REDUC4 can be run, 1f£ desired, to plot a data run on a 
semi-log scale with a user defined D to compare with the 


oe 


meeversal curve. 
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iy vl Late ALK HEATER 


A schematic of the vitiated air heater system at the com- 
bustion laboratory is shown in figure #12. It was designed 
to supply air at any desired mass flow rate within the system 
limitations (approximately 4 lbm/sec) and at any temperature 
up to approximately 1500°R. The heater output air is routed 
to the input sections of the Solid Fuel Ramjet or the Turbo- 
fee combustor. 

The vitiated air heater is fueled by gaseous ethylene and 
bottled oxygen replaces that consumed in the combustion pro- 
cess. Large high pressure air tanks act as the primary air 
supply and the mass flow rate is determined using a sonically 
choked converging nozzle. 

The data acquisition system was designed to monitor and 
control the operation of the air heater in order to provide 
the desired air flow rate and temperature and to provide 
built-in system safety. 

Two Hastings precision flowmeters and two Hastings flow 
controllers were installed to accurately meter and control 
the oxygen and ethylene flow rates to achieve the desired 
fuel-air ratio and oxygen make-up flow rate. Figure #13 
Shows the interconnection of the flowmeters and controllers. 

A description of the system operation is presented below. 


The 'HEATER' program disc is loaded into the HP-85 data 
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Begure 13, Hlow Control Block Diagram 
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Meenisition system. The “SET UP" option is selected to ini- 
tialize variables and perform preliminary calculations to 
determine initial manual control settings. Figure 14 depicts 
the flow chart for the preliminary portion of the 'HEATER' 
Eeogram. 

The operator inputs the desired mass flow rate of air in 
lbs mass/sec. Other inputs are high pressure air tank pres- 
sure (0-3500 PSIG), the estimated cold air temperature (°F), 
desired hot air temperature (°R), and ambient pressure (inches 
Mg). 

The required fuel/air ratio is then computed for ethylene 
and air from a regression of data obtained from the "PEPCODE" 


chemical equilibrium/adiabatic combustion Program then ao |. 


Tyot 7 2°° _ "pus 


F = : 
78000 MatR 





Once F is obtained, the ethylene and oxygen flow rates 


are given by: 


(4 mass/sec) 


For the complete combustion of ethylene with oxygen in air: 


meee Onwt -2H.0 


CoH, 


Zz 
Thus, the desired mass flow ratio of oxygen to ethylene is 


aes = eae 9) Ma 


50 





KEY LABEL 
SELECT OPTION 
HELP SE aur 
DOCUMENTATION NATE 





Compute 
Compute ff 
Compute 
Compute 
Compute 






"FLOW" 





1) Air valve--ON 

2) Measure Ea and aay 

3) Compute mata 

4) Compute APT, to get Marz desired 
>) Set Vo9 

6) Compute mo, and men, 
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The volume flow rates of oxygen and ethylene in SCFM are: 


%, = 734.89 - Mo, 
Qc, = 339.87" - MoH, 
thus 2%, = 300 H, 
where Q is the volume flow rate. The oxygen and ethylene 


flowmeters produce 0-5 Volt DC signals linear with volume 
flow rate (0-20 SCFM for oxygen and 0-6.5 SCFM for ethylene). 
The required oxygen and ethylene set voltages are computed 


~GOM 
Wen = 1927230- M, (lbm/sec) (millivolts) 
2 


V = 646054 -m 
oH, CoM, 


(lbm/sec) (millivolts) 

0), is generated by the dual channel digital to analog 
(D/A) converter card in the 3497A and sent to the REMOTE 
Permminals of the oxygen flow controller (Fig. 13). This flow 
controller is set to the 'BIAS‘' mode of operation with the 
Bocal Set Point (LSP) control set to Zero. 

[echis position (Fig. 15a) the input stages of the flow 
controller sum the three voltages and generate an error 


Sagal : 
BRROR =e UNEUts— 9 bsP- REMOTE 
This signal drives the flow control until the error signal is 


equal to zero, which occurs when the desired oxygen flow rate 
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Buouee 15; Flow Controller Input Diagram 
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is achieved. With the LSP set to 000, the oxygen set voltage 
from the D/A converter acts aS a remote set point to control 
the desired flow rate. The upper section of figure 13 shows 
the interconnections of the second flow controller, the ethyl- 
ene flowmeter, and the oxygen flowmeter in the RATIO mode of 
Operation. In this mode, the flow controller acts to provide 
the proper ratio of ethylene to oxygen under all conditions 
Buecontrolling the ethylene volume flow rate to maintain a 
Meoeratio. Thus, the computer controls the oxygen flow rate 
which in turn controls the ethylene flow rate. 

Mmequre i565 shows the flow controller input circuit in the 
1-5 RATIO mode. In this case, the REMOTE (OXYGEN) signal is 
Meme tO the product of the Local Set Point (LSP) and the 


INPUT (ETHYLENE) signal. 
LSP - INPUT = REMOTE 


However, the input (ethylene) and remote (oxygen) flow 
Capacities are different. Dividing by the full scale flow 


capacities in SCFM normalizes the equation. 


( INPUT ) _ ( REMOTE ) 


LSP “Ss 20 


REMOTE = Oxygen volume flow rate 
INPUT = Ethylene volume flow rate 


mhaen, 


INPUT _ ey 1 a 
REMOTE 20 /\LSP 3 


moms, LSP = 0.975. 
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For example, if the desired oxygen flow rate is 3.5 SCFM 


or 0.00476 lbm/sec, then the OXYGEN set voltage is 


a a) VOLES 


= 0.875 Volts 


The computer sends 0.875 volts to the remote terminals of 
the oxygen flow controller. The RATIO controller responds by 
driving the flow control valve until the ethylene flowmeter 


reads: 


PNP = £5? = "REMOTE 


sO ears 
INPUT = 0.075 7 0.8974 VOLTS 


uw 





ui 


This corresponds to an ethylene volume flow rate of 


, _ 6.5 SCFM 
Qa y= 0.8974 Volts rEUe LCS: 
274 
Q = 1.166 SCFM 
CoH, 


which is 1/3 of the oxygen volume flow rate. 

This completes the preliminary calculations and the sys- 
tem is now turned on to provide air flow without any heat 
eometon. Thus, the operator selects the 'FLOW' option (Fig. 
14) and the computer measures ambient pressure then opens the 
air valve, allowing air to flow to the heater while other 
components are maintained in the off position. High pressure 
air passes through a choked flow nozzle, where Pry and TT, 
are measured. Based on these measurements the actual mass 


flow rate is computed from: 











"AIR ~ = 
i ly 
Jo 
where 
Ph = 
Y = 1.4 AIR 
x2 
at = TS (in?) 
a* = .235 (in) typical 
Ch = 0.97, discharge coef. 
7 feo: 
1.) lbmeR- 
Try s=emeasubed temperature, °R 


Pry = measured pressure, PSIA 


LDR Et 


g aes 2. 2 5 
lbm sec 

The actual mass flow rate of air 
dated oxygen set voltage is computed 
pressure correction APD, is computed 
mass flow rate. 


The operator may then accept the 


proceed to another program segment. 


(lbm/sec) 


is displayed and an up- 
and set. The required 


to meet the desired 


EVOWw COnNalt1Ons and 


If he accepts these 


flow conditions, the actual mass flow rate of air is set 


equal to the desired mass flow rate. 


meen existing flow conditions, Ete 
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Le he -1sSmeee Satistled 


is adjusted manually bv 











applying the recommended APr,.- Tie “FLOW: option is -auco-— 
matically executed again until the desired accuracy is 
attained. 

Once the desired mass flow rate of air is obtained, the 
'START' option is executed. Figure 16 shows a flow chart of 
this program segment. This routine first closes all valves 
and measures Pry which represents ambient pressure under no 
flow conditions. The proper oxygen set voltage is set ac- 
cording to the latest flow calculations. The computer then 
turns on the high pressure air flow, energizes the igniter 
system, and turns on the ethylene and oxygen flow to the 
heater section. The output hot air temperature is monitored 
by a thermocouple and measured by the computer at 0.5 second 
intervals during the start up sequence, beginning 3 seconds 
after ignition. The igniter is turned off after 5 seconds of 
operation unless secured earlier by manual control. 

During heater start up a typical temperature-time trace 
would be as shown in figure 17. 

Temperature is closely monitored during the start se- 
guence and igniter shutdown is executed: 


1) when TyOT AI reaches 700°R (if not previously secured 
by the 5 second timer) 


2) if Tyor arr falls by 50°R in any 5 second interval then 
HEATER shutdown is called to close the OXYGEN and ETHYL- 
ENE valves, close igniter valves, set O7 set voltage = 9 
and close the high pressure air valve. 


= o 7 = 
When TuoT AIR reaches T DESIRED 100° the Heater Mainten 


ance routine is called, terminating the start sequence. The 


ae 
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Paeure 6. “HEATER e Flow Chart=-=Part Two 
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Heater Maintenance Program is assigned to timer #3 and set to 
automatically execute every 20 seconds, while another project 
program segment is run. 

Figure 18 shows the flow chart for the Heater Maintenance 
Program segment (HTRMNT), which is assigned to Timer #3 to 
interrupt the computer at 20 second intervals to control the 
heater output temperature. On interruption, the computer 
measures the output air temperature and compares it to the 
desired alr temperature. If the difference is greater than 
200°R the heater shutdown program is called and secures the 
heater entirely. If the difference is less than 100°R the 
computer records the air temperature and returns to whatever 
routine was operating prior to the interrupt. However, if 
the difference between actual temperature and desired tempera- 
ture is over 100° and less than 200°R the sign of the error 
is determined and the fuel air ratio is adjusted to compen- 
sate for the temperature deviation. Oxygen and ethylene flow 
rates are updated and a new oxygen set voltage is sent to the 
oxygen flow controller to effect a change in the output tem- 
perature. The output temperature is then recorded and a re- 


meen from interrupt is executed. 
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megine IG. HEATER Flow Chart=-Part Three 
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Vee oC lle ar UE RAMG ET 


The Solid Fuel Ramjet (SFRJ) test facility at the combus- 
tion laboratory has been constructed to facilitate the mea- 
surement of performance as affected by design and operating 
variables. Figure 19 presents a schematic of the SFRJ motor 
which is mounted on a thrust stand. Figure 20 presents a 
schematic of the air supply system. Air flow is directed to 
the head-end and sometimes also to the aft mixing chamber 
(bypass-air) of the combustor. In general, a test sequence 
is as follows: air flow rate(s) are turned on, the igniter 
is actuated, a small amount of ethylene is mixed with the 
head-end air flow, the motor ignites, the igniter and ethyl- 
ene flows are terminated, the motor runs for 10-45 seconds 
duration, air flow is terminated and nitrogen gas is flowed 
through the fuel grain for five seconds to ensure extinguish- 
ment of the flame. During the test it is important to mea- 
Sure air flow rates, inlet air temperature, head-end, and 
mixing chamber pressures. Postfire data collected include 
fuel grain weight loss which is needed to calculate the 
average fuel flow rate and average fuel regression rate. 

Temperature rise combustion efficiency for the test is 
determined from 


_ Ttexp 7 Ttair 


" — = 
ae aoe eae 
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where Tt.:+ is the measured inlet air stagnation temperature, 


Teen is the theoretical combustion temperature determined 


from adiabatic, equilibrium combustion calculations using 
the measured fuel-air ratio, inlet air temperature and aver- 


age chamber pressure. Tt is the experimentally determined 


exp 


combustion temperature. The latter 1s calculated from con- 


servation of mass, using gaS properties determined from the 


Sea brium calculations. MThus, 
a 2 vy+1 
T Oepee seh oka a ie 
Cexp my, R =o) 
where Mn = Mir + Mevel 
Aly = nozzle throat area 
Po = average chamber stagnation pressure 
Cy = nozzle discharge coefficient 


Average fuel regression rate is calculated based upon 


fuel weight loss during the burn. Thus, 


=~ 4 AW + 2 
Te mLo ce 


and 


where AW = fuel weight loss 


L fuel grain length 





o = fuel density 
di. =eGinal average fuck internal diameter 
th = burn time 

~ = average fuel regression rate 


The average chamber pressure (Po) used to determine Coe 
is obtained from the pressure time trace recorded during the 
test. A typical trace is shown in figure 21. 

The test control/data acquisition/data reduction program 
that performs the above experiment is described in this 
section. 

Figure 22 depicts the flow chart for the Solid Fuel Ram- 
feces SERJ' option, which 1S assigned to special function key 
#3 on the HP-85. This option is executed after having loaded 
and run the 'HEATER' program (figure 14), completed the 
'SETUP' option, and executed the 'FLOW' option to set up the 
proper flow rates of air, oxygen, and ethylene. If the 
vitiated air heater is to be used, the 'START' option (fig- 
ure 16) is executed to bring the heater up to stable operat- 
ing temperature. Once the air temperature and total air flow 


rate (m,.) are within desired limits of the required values, 


A 
the heater maintenance 'HTRMNT' (figure 18) is automatically 
set to interrupt the computer at 10 second intervals to 
meek the operation of the heater. 

Once the 'SFRJ' key is pressed, the operation of the 


Solid Fuel Ramjet is initiated. Both air flow control valves 


are initially opened to approximately the 50% position. 
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Initially, heater output air is routed via primary and bypass 
air lines to air dump valves, which are held in the dump po- 
sition, venting heater air to atmosphere. The SFRJ igniter 
and the ethylene valves are also maintained off at this time. 

The primary valve position is then adjusted by the com- 
puter until the primary air flow rate fe US Within Strot 
the desired value. The computer programs the 6942A, TIMER/ 
PACER, A/D converter, and memory cards to read the SFRJ motor 
chamber pressure (Pa) transducer output voltage continuously 
for 60 seconds after being externally triggered. The reading 
rate is fixed at approximately 16 readings per second. 

When the 6942 is programmed the SFRJ manual ignition 
switch is enabled by the computer which then monitors this 
switch until it detects closure. Once the switch is closed 
and sensed, the Timer/Pacer card is triggered and the com- 
puter timer is set to zero. Heater air is then directed 
into the SFRJ motor, the igniter is turned on, and the CoH, 
1s turned on. 


The computer waits a preset time in seconds (T ) before 


IGN 
@i@icting off the igniter circuit and the CH, valve. During 
the start sequence the 3456A digital voltmeter is set to 
monitor chamber pressure (Po). The chamber pressure 1s mea-~ 
Sured and compared to the desired steady state value. If 
Maes pressure exceeds .75 of the desired value by the end of 


ignition the computer proceeds to a measurement sequence of 


Several pressures, temperatures, and flow rates until the 
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desired burn time in seconds expires. However, if the pres- 

sure is less than .75 of the desired value the program di- 

rects the user back to attempt another ignition sequence. 
The automatic measurement sequence proceeds from the 


point where the pressure exceeds .75 P as (refer 


C desired 


me taigure 21): 


a) Pry, TT, are measured 











Pp, AC, F 
b) m . ae ee 
air Q 
where Cy = .97 
= na-/4 
Gd Sa3s2 35 
ai 
y=1.4 
aoe 
QO = 
Je 
e) Pros Tp, measured at primary air orifice 
meee. 12 82 Sp Fa 
) — = 
: R 
Dae 
where Q = 
Se 
e) Ving is set to obtain desired Mm. 
f) M™sypass ~ ™a ~ ™pRimary 
Maypass = 0 if no bypass selected 
g) Thor AIR measured 


ak 





my «|OCOT measured 


BP 
i) SFRJ head pressure and temperature measured Piss Ty 
me. “oa “sp +p 
— 
air 


k) Values of time, mair, T, and Tyorp are stored in an 
array for later use in reduction of data. 


The above measurements proceed for a total of 100 read- 
ings at a rate of approximately 1.5 per second. 
Both the SFRJ and the heater are shut down automatically 


when timer #1 expires at (t,) seconds after successful 


B 
iemLcion. 

Once the measurement sequence completes and the shutdown 
sequence is completed the computer then waits until the multi- 
programmer memory board signals completion of its storage 
operation. 

During the operation of the SFRJ option, special function 
key #4 1s assigned to the SFRJ shutdown routine. This pro- 
gram is automatically executed th seconds after lgnition, 
however, it can be invoked at any time by the operator and 
when called performs the following steps: 

a) Air valve to dump position 
b) Nitrogen purge on for 5 seconds 
c) Heater shutdown sequence. 
After all timed events have completed the computer auto- 


matically loads and executes '‘'SFRJ2', which is the data 


WZ 





reduction program for the Solid Fuel Ramjet experiment. Fig- 
ure 23 shows the flowchart of the data reduction program. 

A large quantity of measured data is automatically passed 
from the Heater program to the reduction program. The com- 
puter then prompts the user for other required input data, 
such as: 

a) fuel type HTPB/PMM 

b) final grain weight (grams) 

c) initial grain weight (grams) 

d) initial grain internal diameter (cm) 

e) final aft-end internal diameter (cm) 

f) grain length (cm) 

g) area of the exhaust nozzle throat (em*) 

With this data entered and stored, the computer then reads 
in the chamber pressure (Pa) array, which is a 1000 element 
vector of transducer voltages. These voltages are subsequent- 


ly converted to pressure (PSIA) by the formula 


Po (tT) = K, (v(t) oy a so 
where K, = 79.434 
V(I) = transducer voltage 
V@ = initial transducer voltage at zero psig 
B9 = atmospheric pressure 


The computer then locates the maximum pressure reading 
and its location in the array. The computer then searches 


een forward and backward in time from the location of the 
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Figure 23. Solid Fuel Ramjet Reducticn Flow Chart 
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maximum pressure until it finds the time where pressure has 
dropped to 0.2 times the maximum pressure and records these 


locations as ee: and Cece The effective final burn 


time 1s tap = tp - ta. The computer then performs a numerical 


integration of the pressure time trace from t, to t, and 


divides the result by t,, to give Pa, the average pressure. 


- 
: Po(t)dt 
S 


Pp. = 
E on 


Now that t. and te are known, the computer can locate 


these times in the time column of the data array (100x3 


matrix) and note the row numbers of I. and I... The average 


> 
alr temperature ee and flow rate Moir are then computed 
between to and te by: 
ae 
A 
_ 7 1=1, 
air (I,-I,) 
Tp 
a 
ene 7 1L=1, 
av (I-I,) 


Data reduction can now proceed as follows: 


AW, = Initial W, - Final WwW. 
: 4 AW, 5 
B  \ioeouel a 
P 
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where pp = Fuel propellant density (gms/cm?) 


L length . (cm) 


d. initial diameter (cm). 


ae 


iI 


The average burn rate is then 


. 
Fh 
bv 


Thus, fuel/air and equivalence ratios are 


il 
aU 
ne tr 


and 


£ 

oS (a 

SLOLren 
Ge@enewith T, 6, and the type of fuel burned, the theoretical 
values of Thin? Ras and Ye are computed from a polynomial 
regression obtained from the PEPCODE [Ref. 5] computer pro- 
gram at the Naval Postgraduate School computer center. The 
effect of Be on Tt,, is negligible for nominal variations 
(40-60psia). However, Be effects can be readily incorporated 


if desired. These polynomial equations are given in Appendix 


i 


v7 





Once the values of Tten? Rar and Yq are known then: 


lL. 
\ - 2 
lo 1 


Yacwl 
op ee) emer mez 
Z R Ye t 1 





C 


and the experimental value of chamber stagnation temperature 


Teexp 1s determined from 
Tene J Pe An “p £2 ‘ 
(Mave iy mp) 
where An = throat area 
Cy = .97 


Finally, the efficiency is computed 


a eee T 
aE — 


Hard copy printouts of all data is then provided by the 


system printer. 
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APPENDIX A 


LS teO pee ADDRESSES 


Cet (interna b} 

Printer (internal) 

2631B Line Printer 

Gh (2 Sy) JS oye elene 

3497A Data Acquisition/Control Unit 
He=eo Control ler 

3456A High Accuracy Digital Voltmeter 
6942A Multiprogrammer 

3437A High Speed Digital Voltmeter 
HG-78 Scanivalve Interface #1 

HG-78 Scanivalve Interface #2 

Disc Drive #0 

Disc Drive #1 


Internal Tape Drive 
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APPENDIX B 


PEPCODE REGRESSION EQUATIONS FOR HTPB AND PMM 


Fuel Type £ = 017 3 2 


SLOLCH Pp = .0332 #m/in 


ate eo Po = 40-60 psia 
STOICH 


moons + A6eao 4 2030- = 8966" = 205.90° 4 1.1925T, 


aan See = a 
Te OmNlon st = = 1484 (10 ")T,” = 1.138T, + .3864T,¢ 
eae 
- 7.67(10°7)T,¢ 
2 3 6 
Onc se 4eso- =| .0244e6- + .005236 
: 4.994(10°>)T, + 4.255(10°7)T.7 ~ 2.36 (10713)F.? 
A A 
+ 7.7585 (10° °)T.6 - 3.1138 (1077) T, 6° : 1.885(1077)T “4 
2 3 6 
eoGee acoee roo 5 716°? + 227396 
- af ae eee 
8.08(10°°)T, + 1.36(107/)T, Seo Oe a, 
+ 1.0208 (10 “)T.¢ : 1.872 (1074) F, 4° Z 2.185(107) 7,76 
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zs es ae 
PMM = Fuel Type f mOICH = 0.12 2p = .0426 #m/in 


Pa = 40-60 psia 


= 1642 - 59026 + 146889 - 747847 + 4.057. 6 


Tr 
TH 
= _ 9 Bi eae) 
- .18T, - 4.687, 97 + 1.287, 9° - 3.259(10°>)T,%¢ 
: iG 2 3 
Yo = 1.34 - .0265¢ - 9.126(107°)T, - 19397 + .13296 
Z 4.205(107-)T.¢ + 7.429(107°)T, 4? Z 7 i GON ag 
+ 2.09(10°7)T,* 9° 


R= 47.1569+9[27.54+9{2.9108-9 (41.242-9[19.545-11.5675°])}| 


3\—q 2) 


i (107°) F, {12403 - 3.2029 (107 i f 
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APPENDIX C 


PROGRAM LISTINGS 


eT 


my! HELP FOR HEATER 

20 OPTION BASE i 

30 COM V¢(100,4) ,M9,NOSTi] ,PD,V9,B9 

40 CRT IS 2 

SO DISP "The HEATER program contains project control an 
d" 

60 DISP " DATA ACQUISITION and DATA REDUCTION programs 
for the VITIATED AIR HEATER. " 

70 DISP “PRESS ’SET UP’ KEY TO MAKE PRELIMINARY AIR FLO 
W CALCULATIONS, " 

80 DISP "WHEN YOU ARE SATISFIED WITH YOUR CALCULATIONS 
YOU CAN FLOW THE AIR SYSTEM TO CHECK, " 

90 DISP 

100 DISP “PRESS ’FLOW’ KEY TO COLD FLOW THE AIR TO MEAS 
URE Pti AND Tti AND COMPUTE m:dot:AIR., * 

140 DISP “IF YOU ARE NOT SATISFIED WITH Pti SETTING OR 
ACTUAL AIR FLOW THEN ’?FLOW’ IS RUN AGAIN." 

120 DISP “AFTER ’FLOW’ IS EXECUTED THE ’START’ OPTION K 
EY IS PRESSED TO START HEATER SEQUENCE. " 

£30 DISP "’START’? BRINGS HEATER UP TO DESIRED TEMPERATU 
Re.” 

440 DISP "WHEN DESIRED TEMP IS ATTAINED THE ’HTRMNT’ HE 
ATER MAINT ROUTINE IS RUN EVERY 20 SECS." 

i150 DISP "THE HEATER IGNITER IS DE-ENERGIZED AFTER 8 SE 
CONDS GF OPERATION OR WHEN 700 DEG [IS * 

460 DISP "REACHED , ALSO IF THE AIR TEMP DROPS BY SO DE 
GIN S SECONDS DURING START UP THE IGN" 

470 DISP " IS SHUT GFF AND THE HEATER IS SECURED." 

180 DISP "WHEN THE AIR TEMPERATURE IS WITHIN 150 DEGREE 
S OF THAT DESIRED, THE HEATER MAINT ” 

490 DISP " ROUTINE IS CALLED. THIS PROGRAM EXECUTES EV 

ERY 20 SECONDS UNTIL ANOTHER OPTION" 

200 DISP “ IS SELECTED, MEASURING AND RECORDING THE AIR 
eMP. “ 

210 DISP " THE SOLID FUEL RAMJET OPTION KEY ’SFRJ’ IS P 
RESSED TO ENERGIZE THE SFRJ WITH OR “ 

220 DISP "WITHOUT THE HEATER IN OPERATION. THIS OPTION 
RUNS FOR A PRESET TIME BEFORE SECURING" 

230 DISP " ’°’SFRJ2’ THE DATA REDUCTION PROGRAM IS AUTOMA 
TICALLY CALLED AND EXECUTED AFTER ’SFRJ’" 

240 CRT IS 14 

250 CHAIN "HEATER " 

260 END 
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40 ! HEATER 

20 OPTION KASE 1 

30 COM V(100,4) ,M9,NO$011,P0,V9,R9 
40 ON KEY# 141,"SET UP" GOSUB 150 

SO ON KEY# 2,"FLOW" GOSUR 790 

60 ON KEY# 3,"SFRJ" GOTO 2271 
yaeON KEY#*7,"°TITC" GOSUB i220 

80 ON KEY# S,"HELP" GOSUB 1240 

90 ON KEY# 6,"“START" GOSUB 1270 
400 ON KEY# 4, "SHUTDOWN" GOSUH 2144)) 
{40 CLEAR @ KEY LABEL 

meo DISP "SELECT OPTION" 

430 GOTO 130 


140 ! Initialize the constants as needed to do the prel 
iminary 
450 ! calculations for setting up the air heater. 


460 DISP "WILL YOU BE RUNNING THE HEATER ONLY? (Y/N)° 
470 INPUT Y$ 

180 IF Y¢$="N" THEN GOTO 2880 

190 |! 

200 OUTPUT 709 ;"D04,4,5,6,7" 

210 LET V9=280 ! STORAGE TANK VOLUME: 

220 G=1.4 ! AIR GAMMA 

230 DISP “ENTER STORAGE TANK GAGE PRESSURE (PSIG)" 

240 INPUT H®8 

250 LET C=.97 ! DISCHARGE COEF. 

260 LET Gi=32.2 ! GRAVITY 

270 LET D=1.5 ! INCHES 

280 Di=.235 ! Dstar choked venturi 

2970 LET R=53.3 ! GAS CONSTANT FOR AIR 
300 CLEAR 

310 DISP "ENTER DESIRED MASS FLOW RATE (#/SEC) M9" 

320 INPUT M9 

330 PRINT “ FOR DESTRED MASS FLOW RATE OF "3M9;" #/SEC. 


340 ! DISP "ENTER NOZZLE DIAMETER (Di=COLD AIR=Dx) "“ 
350 ! INPUT Di 

360 DISP "ENTER COLD AIR TEMPERATURE (DEG=F) Ti" 

370 INPUT Ti 

380 LET Ti=T1+460 

390 DISP "ENTER DESIRED HOT AIR TEMPERATURE (DEG=R) TS" 


400 INPUT T8 

410 DISP “ENTER BAROMETRIC PRESSURE SETTING (INCHES HG) 
(29.92=STANDARD SEA LEVEL)" 

420 INPUT B9 

430 B9=,4A9L3XKKY 


440 ! 9 WILL BE USED TO CONVERT ABS PRESSURES TO GAGE" 


450 ! NOW COMPUTE VISCOSITY FOR LATER USE" 
460 LET U=.0000000229*T8*1,S/(T8+198. 6) 
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470 
480 
499 
Soo 
ts 
910 
$20 
530 
$40 
550 
560 
570 
580 
590 


' 
! COMPUTE FUEL AIR RATIO FROM PEPCODE EXCERT 
LET F=(T8-550)/78000 
! This F was for air at approx. 70 degrees farenhei 


! compute ethylene flow rate 

LET M8=FXM9 ! MDOT FUEL 

M7=96/28*M8 ! OXYGEN FLOW RATE 

V=4183720KMN7 ! OXYGEN SET VOLTAGE. 

VS=VALSE(V) 

OUTPUT 709 3"A0D3Z,0,"AV4 

CLEAR 

DISP USING $90 3 N8;826.06*KM3 

IMAGE “SET ETHYLENE FLOW RATE=",1D.9D," #/SEC" ,SX,2 


bee,” SCH" 


600 DISP USING 640 ; M73;722,.S58xM7 

640 IMAGE "SET OXYGEN FLOW RATE=",4D.SD," #/SEC",SX,2D. 
3D," SCFM" 

620 PRINT "OXYGEN VOLTAGE= "3V;" MILLIVOLTS" 


630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
7390 
790 
800 
810 
820 
330 


LET A=PI/4xD1i*2 

LET Fi=SQR(CK(2/(G44))*C(G44)/(G-4£))) 
PRINT “F(GAMMA)="3F4 

LET Q9=SQR(TLXR/G1) 

LET P6=M9*Q9/(CKAKF1) ! PSIA 

LET PS=P6-B9 ! PSIG 

DISP USING 700 ; P&S 

IMAGE “SET PT TO ",6D," PSIG" 

LET TS=(BO-PSK2)KV9K144/ (RET LXM9) 
DISP USING 730 ; TS 

IMAGE "RUN TIME IS ",SD.D," SECONDS’ 
COPY 

CLEAR 

DISP “SELECT OPTION" 

KEY LABEL 

RETURN 

1 FLOW OPTION 

QUTPUT 709 ;"DO04,4,5,6,7" ! AIR VALVE CLOSED 
CLEAR 

WAIT 2000 

! THIS SEGMENT FLOWS THE SYSTEM COLD TO DETERMINE A 


CTUAL MASS FLOW RATE 


840 
650 
860 
870 
880 
890 
9090 
910 
me 


OUTPUT 709 3;"ACS" ! MEASURE PT4 
OUTPUT 722 ;,"HSMODe2SWEZOSOILA  ASTISOFLRITSQX1i® 
GOSUB 2230 


VO=V 

PRINT "VO=",V90 

OUTPUT 709 3;"DC4,6" ! AIR VALVE OPENS 
WAIT 10900 


OUTPUT 722 3; "HSMOO2SW2Z0SOiL1  ASTISOFLRATIQXi" 
GOSUB 2230 
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930 
740 
750 
950 
970 
980 
990 


PRINT "V=",V 
P=79,.434K(V-V0) | 
PRINT "PTi=",P 
OUTPUT 709 ;"AC6" ! MEASURE TT1 

OUTPUT 722 ;"HSW2SM002Z0SO1L1. 1STISOF IRI T1QX4" 
GOSUB 2230 

T=41575kV+536.6 ! VOLTS TO DEGREES R 


Volo tO PSIG 


1000 
10190 
1020 
1030 
1040 
1050 
1060 
1070 
1080 


1090 
1100 
1110 
4120 


(AIR) OF 


1130 


OUTPUT 709 ;"DO4,6" ! CLOSE AIR VALVE: 

1 CONVERSIONS ON T AND P NEED TO BE DONE. 
Q9=SQR(RKT/G1) 

M=CKAXKF1X(P+H9)/Q9 ! ACTUAL MA FLOW 

DISP “ACTUAL MASS FLOW RATEE IS = ="3M;" 
MB=FKM @ M7=96/28%MB 

V=183720%M7 ! OXYGEN SET VOLTAGE 

QUTPUT 709 3"AQ3,0,"AVALS(V) 
P7=M9KQ9/(FL*CKA) |! PSIA 


¥m/sec," 


P=P7~-P-B9 ! PRESSURE 
PRINT “ DELTA PTi=",P 
DISP USING 4120 ; P,M9 
IMAGE “CHANGE Pti KY ",4D," 
",D.5D," #mass/sec," 
"ARE YOU SATISFIED WITH THE MASS FLOW RATE A 


INCREMENT 


PSIG TO ACHEIVE A Mdot 


DISP 


ND DOME PRESSURE SETTING? (CY/N)" 


1140 
L150 
1160 
m7 0 
1180 
1190 
1200 
aye 1 0) 
aa 
1250 
1240 
1250 
1260 
1279 
1280 
1290 
i300 
1310 
i320 
1.330 
1340 
4350 
1360 
a7 | 
1380 
1390 


INPUT Y$ 
IF Y$="Y" THEN GOTO 4180 

Ti=T ! P6=P7! SAVE LAST PT4 AND TT4. 
GOTO 790 

DISP "SELECT OPTION® 

M9=M ! ACTUAL=DESIRED MDOT 

KEY LABEL 

RETURN 

HEEP 200,100 

RETURN 
| HELP 
CHAIN 
RETURN 
! HEATER START UP ROUTINE 
OUTPUT 709 ;"DC4,6" ! AIR ON 
CLEAR 
QUTPUT 709 
OUTPUT 722 
GOSUB 2230 
P=79,434K(V-V0) | 
OUTPUT 709 
OUTPUT 722 
GOSUB 2230 
T=4457S5KV+536.6 |! 
Q9=SQR(R¥T/G4) 
M=CKAKF4K(P+4E9)/Q9 | 


"HELP “ 


pemes6!l OP th 
;"HSWESMO0eCZOSOALA .ASTISOFARATAIQXi" 


VOLTS TO PSIG 

,;"AC6" |! Tri 

; "HSW2SM002Z0S0i1L4.1STISOFAR4TIQXi ” 
T COLD 


ACTUAL MDOT AIR 
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4400 
14410 
1420 
14390 
£440 
1450 
1460 
1470 
1480 
1490 
1500 
45190 
{520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
16090 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1699 
1700 
4710 
4720 
4730 
4740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
4830 
1840 
1850 
1860 
4870 
1880 
1890 


MB=FaM @ M7=96/28KMB 
YV=4{B3720eMN7 ! Oe SET VOLTS 


QUTPUT 709 ;"A03,0, "AVAL$(V) 

OUTPUT 709 ;"AC7" ! T hot 

OUTPUT 709 ;"DC4,7" ! IGNITER ON 

OUTPUT 709 ;"DC4,4,5" ! O02 AC2H4 ON 

ON TIMER# 41,8000 GOSUB 1670 

T3=T ! T LAST 

WAIT 5000 

QUTPUT 722 ;"HSW2SM002Z0SO1L41.4STISOFiRiT1OQXi" 
GOSUB 2230 


T2=1002. 34U+457.53 |! 
T1=T2-T3+S0 

IF T2)700 THEN GOSUB 1720 

IF Ti¢0 THEN 1760 

PRINT "Tihot=",T2 

1 CHECK TO SEE IF HEATER IS WITHIN DESIRED LIMITS 
IF T2+150-T8>0 THEN 1840 

1 HEATER IS NOT YET NEAR OPERATING RANGE 

OUTPUT 709 3"AC7" 

T3=T2 ! LAST T HOT 

WAIT 4000 
OUTPUT 722 
GOSUE 2230 
T2=1002.34V+457.53 |! 
TisT2-T3+S0 

GOTO 1540 ! READ T HOT AGAIN 
1 SHUT OFF IGNITER 

OUTPUT 709 ;"D04,7" 

PRINT “IGNITER SHUT DOQWN-S 
OFF TIMER® 1 

RETURN 

OUTPUT 709 ;"D04,7" 

OFF TIMER® 1 


Thot 


,"HSWeESMODeAZOSOILA A1STISOFARATIQXi” 


Thot 


SES.” 


DISP “ T:HOT > 700 R IGNITER OFF" 

RETURN 

! TEMP GRAD IS NEG TIME FOR SHUT DOWN 
QUTPUT 709 ;"D04,7" 

OUTPUT 709 ;"D04,4,5" 

OUTPUT 709 ;"A03,0,0" 

OUTPUT 709 ;"DO4,6" ! AIR OFF 

DISP " HEATER SHUTDOWN-NEG TEMP GRADIENT* 


OFF TIMER# 14 

STOP 

! HEATER IS NEAR OPERATING POINT 
OFF TIMER# 14 

mi weeitey £74 CLOSE” 

ON TIMER# 3,10000 GOSUH 1949 

DISP "ENTER HEATER MAINT ROUTINE" 
CLEAR 
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7 0 0 
19140 
14920 
1730 
1940 
1950 
i960 


4970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2iid 
2120 
21350 
2140 
2i50 
2150 
2i70 
2180 
2190 
2200 
2210 


2220 


22350 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 


2330 


DISP "HEATER RUNNING" 

WAIT S00 

KEY LABEL 

GOTO i187 

' HEATER MAINTENANCE 

QUTPUT 709 ;"AC7" 

OUTPUT 722 3; "HSMOO2CSWAZOSOILIFLO. OUISTISOFIRITIQXi”" 


GOSUB 2230 ! TAKE READING 
T2=1002.3*V+4S57,.53 ! Thot 
PRINT “T:ihot:AIR=";T2 
T4=TB-T2e ! DELTA T:ihot 

IF ABS(T4)>200 THEN GOTO 2140 
IF ABS(T4)¢400 THEN 1940 

IF T4¢0 THEN GOTO 2060 

DISP “" NOT HOT ENOUGH " 

GOTO 2080 

DISP " TOO HOT " 

GOTO 2080 

F4=7T4/78000 

F=F+F4 ! NEW FUEL AIR RATIO 
MB=FXKXMO @ M7=96/28KMS 
V=183720N7 
OUTPUT 709 
RETURN 

! HEATER SHUTDOWN 

OUTPUT 709 ;"A03,0,0° 

OUTPUT 709 ;"D04,7,4,5" 

PRINT “HEATER SHUTDOWN" 

PRINT " LAST T:hot="3T2 

WAIT S000 

OUTPUT 709 3;"D04,4,5,46,7” 

OFF TIMER 3 

RETURN 

! SRQ SERVICE ROUTINE 

STATUS 7,1 3 Z9 

ENTER 722 3; V 

ENABLE INTR 7;8 @ CLEAR 722 @ RETURN 
! SOLID FUEL RAMJET OPTION 


5"A03,0, "AVAL$(V) 


OUTPUT 709 3"DO4,8,9,49,44,412,43" ! AIR DUMP 
OFF KEY# 4 @ ON KEY# 4, "SHUTDOWN" GOSUB 3050 
IF NO$="N" THEN 2330 

OUTPUT 709 3;"DC4,9,43" 


DIM 00$[90] 
OO$="UF,3.41,0,3.2,0,3.3,0T,SF,2,3,1,.004,12T,MI,2, 


{000T,WE,3,999T,AC,3T,WF,4.2,4T" 


2340 
2350 
2360 


OUTPUT 723 300% 
OUTPUT 723 ;"0P,4,.03ST” 
CLEAR 722 @ OUTPUT 722 ;"HSMO02SW2Z0SO1LiFLO.O4STI 


SOFLR1T3Q" 
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Bove OUTPUT 709 3;"DE4,14" &@ DISP "READY FOR SFRJ IGNITI 
ON” 

2380 OUTPUT 709 ;"ACi3" | Pchamber 

2390 OUTPUT 722 3;"XKi" ® GOSUB 2230 

2400 V9=V '| INITIAL Pc READING 

2440 OUTPUT 709 ;"ACi4" ! IGNITION 

2420 OUTPUT 722 ;"Xi" @ GOSUBR 2230 

2430 IF ABS(V)¢(3 THEN 2420 

2440 ! BEGIN SFRJ START UP 

2450 SETTIME 0,0 

2460 OUTPUT 709 ;"DC4,8,11,12" 

2470 WAIT 1000xT6 

2480 OUTPUT 709 ;")04,12,11" 

2490 OUTPUT 709 ;"ACi3" ! Pchamber: 

2500 OUTPUT 722 3"Xi" @ GOSUB 2230 

pero ITF V),.75x*(P0/151.5) THEN 2540 

2520 CLEAR 

2530 DISP "NO IGNITION TRY AGAIN" ®@ GOTO 2420 
2540 ! REGIN MEASUREMENT SEQ, 


2550 ON TIMER® 4£,76%1000 GOSUER 3050 

2560 OUTPUT 722 ; "HSMNO02SW2ZOSOILAFLO .OLSTISOFIR1ATIQ" 
2570 FOR J=i TO 100 

2580 OUTPUT 709 ;"ACS”" @ OUTPUT 722 ;"Xi" @ GOSUB 22350 
2590 P=79,434kK(V-V0) @ OUTPUT 709 3;"ACS" 8 OUTPUT 7e2e ; 
"Ki*® @ GOSUB 2230 

2600 T=41575*V4+536.6 @ VII, 1LI=CKAKFALX(CP+B9) /SQRCRKT/G1i) 


2610 OUTPUT 709 ;"ACi2" @ OUTPUT 722 3; "Ki" @ GOSUB 22st 


2620 T3=415754V+536.6 @ OUTPUT 709 ;"AC7" @ OUTPUT 722 
;"X4i" @ GOSUB 2230 

2630 V(J,4)=41002.3*V+457.5 @ T2=¥(J,4) ! Tihot 

2640 OUTPUT 709 ;"AC8" @ OUTPUT 722 ;"Xi" @ GOSUB 2230 
2650 P=79.434k(V-V0) @ OUTPUT 709 3"AC9" @ OUTPUT 722 
"X4" @ GOSUB 2230 

2660 T=4157S5%V4+536.6 @ MB=CKAMF 1K (P+B9)/SQR(RKT/G1) 
2670 M4=V(J3,1)-M3 @ YV(J,3)=TIME 

2680 IF NOS="N" THEN M4=0 

2690 OUTPUT 709 ;"AC44" @ OUTPUT 722 ;"Xi" @ GOSUB 2230 


e 
? 


2700 P4=79,434«(V-V0) @ OUTPUT 709 ;"ACiO0”* @ OUTPUT 722! 
;"Xi" @ GOSUB 2230 

2740 T4=41575*KV+S536.6 @ VJ, 2.=(MSKTS44+M4KTS)/UCT,4)> § T 

bar 

e7eud PRINT “J="3;J,"TIMES="3VCI,3) 

2730 PRINT "Thead=",T4 

2740 PRINT “Phead=",P4 

2750 PRINT "T:hot=",V¢J,4) 

2760 PRINT "“m:idot AIR=",V¢CJ,4) 
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2770 PRINT “m:idot PRI=",M3 

2780 PRINT "T;bypass=",T3 

2790 PRINT "T:BAR=",V(J,2) 

2800 NEXT J 

2840 O0=SPOLL(723) ! WALT FOR MEM INTR 

2820 IF 004264 THEN 2810 

2830 DISP "MULTIPROGRAMMER HAS P:chamber data " 

2840 SEND 7 ; UNL MLA TALK 23 SCG 12 

2850 ENTER 723 ; O4@ DISP “ SLOT#=";01 

2860 OUTPUT 723 ;"DC,31T,WF,4.2,0T,0P,4,0T" 

2870 CHAIN "SFRJ2" 

2880 DISP "WILL YOU BE RUNNING THE SOLID FUEL RAMJET? ¢ 
Y/N)" 

2890 INPUT Y$ 

2900 IF Y$="Y" THEN GOTO 2940 

2940 CLEAR 

2920 DISP "TUREBOJET COMBUSTOR OPTION IS NOT AVAILABLE A 
T THIS TIME. SELECT ANOTHER OPTION" 

2930 GOTO 450 

2940 ! SFRJ OPTIONS 

2950 CLEAR 

2960 DISP "ENTER DESIRED CHAMBER PRESSURE (PSIA)" 
2970 INPUT PO 

2980 DISP “ENTER THE DESIRED BURN TIME IN (SECONDS)” 
2990 INPUT TS 

3000 DISP "ENTER THE IGNITION TIME IN (SECONDS) " 

3010 INPUT Té 

3020 DISP “WILL THE BYPASS OPTION BE USED? (Y/N)" 
3030 INPUT NO$ 

3040 GOTO 190 

30S0 ! SFRJ SHUTDOWN 

3060 OUTPUT 709 5"DO4,8,9,414,42,43" 

3070 OUTPUT 709 ;"DC4,10" |! N2@ PURGE ON 

3080 WAIT S000 @ OUTPUT 709 ;"DO04,40" 

3090 GOSUB 2440 ! HEATER SHUTDOWN 

31400 CLEAR @ DISP "“SFRJ& HEATER SHUTDOWN " @ KEY LABEL 
3440 RETURN 
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10 


990 
i100 
1190 
{20 
130 
140 
i50 
160 
170 
i180 
170 
200 
21d 
220 
230 
240 
250 
260 
270 


2890 
299 
500 
310 
320 
330 
Mag 
540 
350 
360 
370 
580 
390 
400 
410 
420 
430 
449 
450 
46) 
470 
480 


1 SFRJ2 DATA REDUCTION 
OPTION EASE 14 
Ome 70100 ,4) M9 NOB 01) ,P0,V9,R9 

SHORT P(1000) 
GCLEAR 

CLEAR 
Gi=32.174 ! gc 

C=.97 ! DISCHARGE COEF. 
DISP "SFRJ DATA REDUCTIONS" 

DISP "ENTER INITIAL GRAIN WEIGHT=(GRAMS)" 
INPUT Wi 

DISP "ENTER FINAL GRAIN WEIGHT=(GRAMS)" 
INPUT We 

DISP “ENTER GRAIN LENGTH=(CM) " 

INPUT L2 

L2=L2/2.54 

DISP "ENTER INITIAL GRAIN DIAMETER=(CM) " 
INPUT DO 

DO=D0/2.54 

DISP "ENTER FINAL AFT DIAMETER=(CM) 
INPUT De 

D2=D2/2.54 

AL=PT*1L.,2*%2/4 ! AREA THROAT 

DISP "ENTER FUEL TYPE=(HTPB/PMM)" 

INPUT Fis 

IF Fis="“HTPEK" OR Fis="PMM" THEN 290 

DISP "“GNLY FUEL TYPES AVAILABLE ARE <HTPH) OR <(PMM> 


GOTO 240 

W3=(Wi-W2)/453.6 ! DELTA WEIGHT IN 3m 

! ENTER Pc DATA 

SEND 7 ; UNL MLA TALK 23 SCG S 

FOR I=4 TO 997 STEP 4 

ENTER 723 USING "#,K,K,K,K" 3 PCI), PC(I+4),P(I+2) , PC 
) 

NEXT I 

ENTER 3 

j 

MAT P=(S) xP 

FOR I=4 TO 1000 

P(I)=15,97KP(1)+.22 1! PSIG 

NEXT I 

Z0=AMAX (P ) 

JO=AMAXROW @ TO0=J0X. 04 

PRINT "MAX Pc="3Z0;" AT TIME=";T0;" SECONDS" 
FOR I=J0 TO 1000 

IF P(I)>.2*Z0 THEN 490 


Se=] 
PRINT "S2=";S2 
GOTO 500 
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490 NEXT I 

S00 FOR I[=J0 TO i STEP i 

S40 IF P¢(I)>,.2*Z0 THEN S50 

S20 Si=lI 

SOuePRINT "Si=";Si1 

S40 GOTO S50 

SS0 NEXT I 

560 S$3=S2-Si ! ELEMENTS INCLUDED 
S70 S4=,06%*Si @ SS=,064S2 

580 S6=,.06*S3 ! TIME TBF 

S90 REDIM P(S3) 

600 S7=SUM(P ) 

640 Pi=S7/S3+B9 ! Pe bar(PSIA) 
5620 REDIM P(1000) 

630 GCLEAR @® LOCATE 20,120,20,80 @ FRAME @ SCALE 0,60,0 
, 100 

e440 LAXES S,100,0,0,2,5,5 

650 FOR I=e TO 1000 STEP e 

660 PLOT I%.05,P(T) 

670 NEXT I 

680 MOVE S4,0 @ PLOT $4,0 @ PLOT S4,P(S1i) 
690 MOVE SS,0 @ PLOT SS,0 @ PLOT SS,P(S2) 
700 MOVE J0*.06,0 @ PLOT J0x*.06,0 @ PLOT J0*,06,P(J0) 
710 COPY 

720 !' FIND M: DOT: AVE AND Thar: AVE 
730 FOR I=i TO i00 

740 IF VCI,3)9<S4 THEN 770 

750 V4=] 

76090 GOTO 7890 

770 NEXT I 

780 FOR IT=V4 TO i100 

790 IF V(1,3)¢4SS THEN 820 

800 VS=I 

840 GOTO 839 

820 NEXT f 

830 V6=0 @ T=0 #6 T3=0 

840 FOR I[=V4 TO VS 

850 V6=V64+1VU(T, 14) !§ M: DOT AVE 

860 T3=T3+V(I,4) ! T:HOT AIR 

870 T=T+VUC1T,2) ! T: BAR AVE 

880 NEXT I 

890 V8B=VS-V4 ! TIME INT 

900 VE&=V6/VB |} M: DOT AVE 

910 T3=T3/V8 ! T:HOT:AIR : AVERAGE 
920 T=T/VB ! Tibar AVE 

930 PRINT "Pc:bar="3;Pi;" PSTA" 
940 PRINT "Tbhf="3;S6;" SECONDS" 
950 PRINT "mm: AVE="3V6 

960 PRINT "T: BARave=";T 

770) PRINT "“Tihotsair=";T3 


oan 





980 PRINT 


990 V7=W3/S6 ! 


4000 
4010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
10970 
1100 
44110 
1120 
1130 
1140 
4150 
1160 
1170 


FUEL FLOW RATE 
X4=V7/V6 ! Fuel/air ratio 
PRINT "f=": X4 
IF Fis="HTPR" 
IF Fis="PMM" 
GOSUB 1420 
GOSUB 1450 
T2=(PLXALKCKF2/(V64VU7))%2 |! Tcexp 

NO=(T2-T)/(Ti-T) !' EFFICIENCY 

PRINT “Tce; THEORY="3; 74 

PRONE "“TC:EXr =.“ Fe 

PRINT "T;: BAR: AVE="3T 

PRINT “THERMAL EFFICIENY ="3N0 

END 

|! HTPB SUBROUTINE 

Be=-.075 ! FrSTOIcCH 

D9=.0332 ! DENSITY #m/in*3 

X=X1i/F9 

fia ay So. Star Sante SkA 2-SI99KA°S=205.9KX%6+1 «1 925k 


THEN GOSUB 4130 
THEN GOSUH 1280 


T+. 00010144KT*%2-, 000000018 4KT*3-1.1358KT KX 
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TI=Tit+,. SB64KTKX*2-. 00007 67%T*2KX 


4490 G=4.40162-.277234%X+. 1 483RX*2-. 02448%X*3+, 00S 2TKX*4 
~,00004994KT+.0000000042S5SKT*2-2, 36E-413T*3 

1200 G=G+, 000077S8SKTX—. 0000341 38%T*X*2—-.00000000188Sk 
T*2%x 

£240 R=S52.966+3,S46*X-9 , SSKX* 245, 74 KX* 34+, 2739KX*6-.0008 
OSKT+.000000436xT*2-5 , 99E—144T*3 

{220 R=R+. 001 0208KTEX-. 0004872KTKX*2+,0000002185%XKT*2 
{230 PRINT "“PHI="5x 
£240 PRINT "Te:th="3T 
1250 PRINT "Re: th=";R 
1260 PRINT "GAMMAc=";G 

{270 RETURN 

1280 ! PMM SUBROUTINE 

1290 F9=,42 ! F:STOICH 

1300 D9=.0426 ! DENSITY #m/in*3 

4340 X=X4/F9 

1320 1T1=1642-5902KX+1 4688KX* 2-74794X* 344, OSKTEX-. L8XT-4 
,6BKTHX* 244, 2BKTHX*4—,00003259KT*24X*2 

1330 G=1.34—, 0265KX—. 193KX* 24+, 4329KX*3—-, 000009426KT-. 00 
00420SKTXX+, 00007429KTKX*2-, 0000272KTKX* 4 

{340 G=G+,.00000000209KT*2*xX*2 

1350 R=47.1569+XK(27.S44XK(2,9408-KK( 44, 242-XK( 19. 545-4 
4.567%X*3)))) 

1360 R=R+.000004%TK(42403-,0032029%T*2) 

{370 PRINT "PHI=";xX 

1380 PRINT "Te:th="5T4 

{390 PRINT “Re:th=":R 

1400 PRINT "GAMMAc="5G 


x 
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1410 
1420 
4430 
1440 
i450 
1460 
1470 
14890 
1490 
{500 
Hoi0 
1520 
ed! 


RETURN 

! FCGAMMA:C) 

FE=SQR (GKGI/R EC A/S (G41) )0*%*0(64+1)/706-1))) 
RETURN 

! D: FINAL: AVERAGE 


DB=SQR (44WS/ (PTXD9XKL2)+D0%2) §' D: FINAL AVE. 


RB=(D8-D0)/(2*S6) !| AVE HURN RATE 
D8=2.54*D6 ! IN CM. 

R8=2,54«eR8 ! IN CM/SEC 

PRINT 

PRINT "d: final: AVE=",D8 

PRINT "AVE: burn RATE="3;R8 

RETURN 


SO A SN we ee ee 
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io ! 
2d! 
30! 
40 ! 
IMENT. 
50 ! 
60 ! 


BILE? 
RKERT HANSEN 


THIS PROGRAM INITIATES THE LIGHT DIFFRACTION EXPER 


IT ALSO PROVIDES HELP WITH THE DOCUMENTATION. 


70 OPTION BASE i 


80 CLEAR 


90 COM Q1$(12)] ,Q@2$(12),T1$l20),DIsliel ,Deasli2l 


{00 DISP 


"ENTER DATE: TIME C(MMDDHHMMSS)D " 


{10 INPUT T1% 


420 OUTPUT 709 
{30 OUTPUT 709 
{40 Q$="DATA"™ ! 
150 CLEAR 


160 F=0 


;"TD"STi4 
, ee ba 
TYPICAL NAME: 


FLAG RESET 


{70 CRT IS 14 


480 DISP 
ON," 
770 DISP 
200 DISP 
210 DISP 
220 DISP 
E50) DISP 
240 DISP 
250 DISP 
he light 
260 DISP 
270 D1ISP 


280 IF F= 


270 DISP 


"ENSURE THAT THE FOLLOWING EQUIPMENT IS TURNED 


sfc nro Pee equrIr. * 

.VISICORDER” 

/.RETICON ARRAY, " 

»LASER “ 

F S.BATTERY AND IGNITION CONTROL CIRCUITS." 
6,E0W PASS —& TLITER™ 


"The following programs are asssociated with ¢t 
diffraction experiment,” 


= 
db Wre 


: PRESS [CONT] WHEN READY." 
0 THEN PAUSE 
*" 4. ACQDTA- acquires data from both diode arra 


¥yS and stores it on the disc." 


300 DISP 
310 DISP 


" 2@.,CMBINE-combines raw data arrays into TWO i 


024x2 arrays and datafiles." 


320 DISP 
330 DISP 


“ S3.PLTDTA-plots the raw data from the datafil 


e as VOLTAGE VS. DIODE NUMBER," 


340 DISP 
550 DISP 


360 IF F= 


370 DISP 


. PRESS ({CONT] WHEN READY. * 


0 THEN PAUSE 
" 4. REDUCE-performs the reduction math and plo 


ts NORMALIZED INTENSITY VS. THETA. " 


380 DISP 
390 DISP 


" S,REDUC2-plots the ideal curve from disc. 


Ca 


libration data can also be plotted." 


400 DISP 
410 DISP 


i“ PRESS (CCONT] WHEN READY.” 


$20 IF F=0 THEN PAUSE 
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meu DLSP =~ 6.REDUC3-Computes D32 from the user specifie 
d datafileand calcs. theta bar." 

440 DISP 

450 DISP " 7,CALIB- If a run is to be a calibration cur 
ve this” 


460 DISP " program performs a curve fit to get a smooth 
curve for comparison." 
Au DISP 


480 DISP “ 8.REDUC4-Plots INTENSITY VS THETA BAR" 

490 DISP “using D32 enterred by user in format like RED 

ine 2 

S500 IF F=0 THEN PAUSE 

Si0 IF F=2e THEN GOTO S8q 

S20 DISP "® DO YOU WANT A HARD COPY OF THE INSTRUCTIONS? 
mi 7N)* 

S30 INPUT A$ 

540 IF AsS="N" THEN GOTO S80 

S50 F=e ! FLAG SET 

S560 CRT IS 2 

$70 GOTO i180 ! RETURNS TO PRINT 

580 CRT IS 14 

590 CLEAR 

600 CHAIN “ACQDTA" 

610 END 
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400 ! ACQDTA 

140 ! BERT HANSEN 

420 OPTION BASE 1 

130 OUTPUT 709 "TOSOO" |! TIME SIGNAL TO VISICORDER 

440 COM Q1$(12],Q2$(421,74$(201 ,D1$(12] ,D2sli2] 

{50 SHORT A( 4032) ,B(4032) 

460 OUTPUT 709 ;"TD" 

{70 ENTER 709 ; Tis 

180 PRINT "DATE: TIME IS ";Ti4 

190 OUTPUT 709 “AR,DO4,0,2,3" 

200 OUTPUT 723 ;"CC,2,3,12,13T,WF,3.4,0,3.2,0,3.3,0T, WE 
,43.4,0,13.2,0,43.3,0T" 

240 CLEAR 

220 CLEAR 722 

230 DISP "ENTER THE FILENAME OF THE DATA FILES TO HE CR 
EATED. (RAWL,RAW2) (MAX 142 CHAR.)" 

240 DISP "DATA SHOULD BE STORED ON DRIVE #4 BY TYPING 
(;D704) AFTER FILENAME. " 

250 DISP “ BE SURE THE FILES NAMED DO NOT ALREADY EXIST 
ON DISC” 

260 INPUT Di$,D2$ 

270 CLEAR 

280 MS=1024 

290 CREATE Di$,1024, 14% 

300 CREATE D2$,1024,46 

340 PRINT " DATA WILL BE STORED ON DISC WITH FILE NAME 
AS FOLLOWS. " 

320 PRINT "NO PARTICLES-" ,Di$ 

330 PRINT "PARTICLES-“,D2$ 

340 DISP "ENTER THE THRESHOLD PRESSURE TO TRIGGER DEVIC 

ES (PSI)" 

350 INPUT S8 ! DESIRED PRESSURE 

360 DISP "ENTER TIME DELAY IN SECONDS FROM THRESHOLD PR 

ESSURES (SECS.)" 

370 INPUT TS 

380 T8=(T8-.25)*1000 

390 V0=S8/154.5 

400 | 

440 CLEAR 

420 DISP “ LASER ON ,IGNITION OFF" 

430 D3$=D1is 

440 GOSUB 840 ! CALL MULTIPROGRAMMER 

450 GOSUH 14180 ! STORE DATA 

460 OUTPUT 709 "ACO" |! CONNECT 3456 TO IGNITION FIRING 

LINE 

470 LOCAL LOCKOUT 7 

480 R7=225 

490 OUTPUT 722 "HSMO02SW2SO1L1SOF4R4T3QXi4 " 

S00 GOSUR 4260 ! ENTER READING 

S40 RB=V 
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$20 TF RB¢CRZ THEN GOTO 570 

S30 ! CONTINUITY CHECK 

S40 DISP "CONTINUITY CHECKS BAD RECHECK BEFORE PROCEEDIL 
NG, THEN PRESS(CONTI" 

S50 PAUSE 

560 GOTO 4990 

S70 CLEAR 

S80 DISP “CONTINUITY IS OK" 

S90 OUTPUT 709 ;"DC4,2" ! NITROGEN Of 

600 DISP 

640 OUTPUT 709 ;"AC2" ! CONNECT DVM TO FIRE SWITCH 

620 DISP “ ENSURE VISICORDER IS SET UP TO RUN ON PROPER 
SCALE AND LAMP IS ON" 

630 OUTPUT 722 “HSMOOG2SW2ZO0SOLILAFLO.OLSTISOFLRIT3Q" 

640 DISP "STANDING BY FOR IGNITION" 

650 KEEP 

660 OUTPUT 7ee "Xi" @ GOSUB 1260 ! ENTER DATA 


670 RY=ABS(V) @ IF R9C10 THEN GOTO 640 

680 OO=TIME @ OUTPUT 709 “DC4,0" @ ON TIMER# 41,5000 GOS 
UB 1250 

690 OUTPUT 709 "“ACi”" @ CLEAR 

700 OUTPUT 722 ;3;"Xi" @ COSUB 1250 ! ENTER DATA 


740 R9=ABS(V) @ IF R9<VO THEN GOTO 700 

720 WAIT T8 @ O1=TIME 

730 GOSUB 840 ! MULTIPROGRAMMER 

740 WAIT 3000 

750 OUTPUT 709 ;"D04,0" 

760 LOCAL 7 

770 D3$=D2% 

780 GOSUB 4180 ! STORE DATA 

790 PRINT “TIME FIRE TO MULTI CALL=";04-00 

800 PRINT "TIME FIRE TO MEM INTR=";02-00 

840 PRINT "DATA STORED ON DISK WITH FILENAMES (";Di$;*) 
AND (";D2$5")" 

820 CHAIN "CMBINE" 

930 END 

840 ! MULTIPROGRAMMER ROUTINE 

850 OUTPUT 723 ;"SF,2,3,4,.004,42T,SF,12,3, 
860 OUTPUT 723 ;"WE,3.1,0T,WF,3.2,0T,WF,3.3 
OT ,WF,43.2,0T,WF,43.3,0T" 

870 OUTPUT 723 ;"WF,3,1034T,WF,43,40347,MI,2,4032T" 
880 OUTPUT 709 ;"DC4,3" ! MARKS THE VISICORDER TRACE 


,.004,42T* 
OT,WF,43.4, 


4 
a 
> 


890 OUTPUT 723 ;"AC,3T,WF,4.2,47,0P,4,16.5T" 
900 K=SPOLL(723) ! WAIT FOR MEM INTERUPT 

940 IF K<)44 THEN GOTO 900 

920 O2=TIME 

930 QUTPUT 709 ;"DO04,3" 

940 DISP “ DATA TAKEN AND STORED* 
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950 SEND 7 ; UNL MLA TALK 23 SCG 12 

960 ENTER 7 ; K4 

970 DISP "MEMORY BOARD HAS DATA" 

980 DISP " SLOT#=";Ki 

990 OUTPUT 723 ;"DC,3T,WF,4.2,0T" 

1000 ENABLE INTR 738 

{040 SEND 7 ; UNL MLA TALK 23 SCG S 

1020 FOR I=4 TO 1025 STEP 8 

{030 ENTER 723 USING "#,K,K,K,K K. mike Kn AX TAC htd 
A(1+2) ,ACI+3) ,ACI+4) ,ACI#5) ,ACI+6) ,ACT+7) 

{040 NEXT I 

£050 ENTER ? 

1060 DISP “ ARRAY ONE DATA ENTERRED" 

1070 OUTPUT 723 3"M1I,12,41032T" 

{080 SEND 7 ; UNL MLA TALK 23 SCG S 

1090 FOR I=4 TO 1025 STEP 8 

{400 ENTER 723 USING "#,K,K,K,K,K,K,K,K" 3 BCI) ,B(I+4), 
BC I+2),H¢1+3) ,BCI+4) , BCI+S) , B(I+6) , BC I+7) 

4440 NEXT I 

1420 ENTER ? 

1430 DISP " ARRAY TWO DATA ENTERRED" 

1140 MAT A=-A 

1450 MAT B=-B 

4460 OUTPUT 723 ;"CC,2,3,42,13T" 

1470 RETURN 

1480 ! THIS ROUTINE STORES THE RAW DATA ON THE DISK 
1190 ASSIGN# 4 TO D3$ : 

1200 FOR I=7 TO 1027 STEP 4 <7 

{240 PRINT# 4 3 ACI), ACI+4),AC1+2) ACI+3) (IT), H(I+4),B 
(142), B(1+3) 

{220 NEXT I 

1230 ASSIGN# 1 TO x 

1240 RETURN 

£250 OFF TIMER# 4 @ OUTPUT 709 ;"D04,0,2" @ RETURN 
{260 ! SRQ SERVICE 

1270 STATUS 7,4 ; Ad 

{280 ENTER 722 ; V 

{290 ENABLE INTR 738 @ RETURN 
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10 !' CMBINE 

20 !' COMBINES RAW DATA FILES INTO A SINGLE DATA FILE 
30 OPTION BASE 1 

40 COM Q1$(421],Q2$(4121 ,T1$(20) ,Di$(12] ,D2e$(12] 

SO D9$=":D70i" 

60 CLEAR 

FO} 

80 SHORT Y(1024,4) 

90 D3$=Dis 

400 FOR M=i TO 2 

1140 ASSIGN# 2 TO D34 

{20 FOR I=i TO £024 STEP 4 

{30 READ# 2 3 YCI,M),YCI+£,M) ,YCI+2,M) ,YCI+3,M) ,YC(I,M+2 
>, YCI+d ,M+2) ,YCI+2 M+2) ,YCI4+3,M+2) 

440 NEXT I 

4S0 ASSIGN# 2 TO 4 

1460 D3$=D2$ 

470 NEXT M 

180 Q3$=Di$ 

190 FOR M=i TO 3 STEP 2 

200 ASSIGN# i TO Q35 

210 FOR N=i TO £1047 STEP 8” 

220 ASY(N,M) @ B=Y(N+4,M) @ CHY(N+2,M) @ D=Y(N+3,™) 
230 E=Y(N+4,M) @ F=Y(N+S,M) @ G=Y(N4+6,M) @ H=YO(N+7,M) 
240 R=Y(N,M+i) @ S=YC(Nt4 , Mti) @ T=YON4+2,M+i) @ UH=YOCN43, 
M+i) @ V=Y(N+4,M+i4) 

250 WHYCN+S ,M+i) @ X=YCN+G6,M+i) @ Z=YCN+7 ,M+4) 

260 PRINT# 4 3; A,B,C,D,E,F,G,H,R,8,T,U,V,W,X,2 

270 NEXT N 

280 ASSIGN# 1 TO xk < 

290 Q3$=D2$ 

300 NEXT M 

340 DISP " DATA STORED ON DISK ON FILE ("3;Did;") AND (" 
‘Deo: ")" 

320 PRINT 

330 PRINT "DATA FILES HAVE BEEN REARRANGED A FOLLOWS: " 
340 PRINT "PHOTO ARRAY #1 IS "“3;D1i$ 

350 PRINT "PHOTO ARRAY #2 IS “;D2$ 

360 CHAIN "PLTDTA" 

370 END 
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40 ! PLTDTA REVISION 1 
20 ! 


30 ! PLOTTING RAW DATA FROM ARRAY 


40 ! 
S50 ! BERT HANSEN 
60 CLEAR 


70 OPTION BASE i 

$0 COM Q1$(412),Q2$(42] ,T1$(201] ,D1i$(412],D2$(121 

90 OUTPUT 709 ;"TD" 

100 ENTER 709 ; Tis 

140 DISP “ ENSURE PLOTTER IS READY FOR PLOTTING WITH CL 
EAN PAPER AND DESIRED PEN." 

120 DISP “PRESS [CONT] WHEN READY!“ 

{30 HEEP 

{40 PAUSE 

150 DISP " ENTER THE FILE NAME OF THE DATA TO BE PLOTTE 
D.(";Di$;" OR ";D2$;" )* 

160 DISP "DRIVE MUST BE SPECIFIED AS (NAME:D704)" 

{70 INPUT Q3$ 

{80 DISP “ENTER DESIRED PLOTTING INTERVAL OF DIODE ARRA 
Y (2,4,6)" 

£90 INPUT J 

200 CLEAR 


210 DISP "READING DATA FROM DATA FILE ON DISC DRIVE 


220 PLOTTER IS 705 

230 SHORT Y(4024,2),X(1024) 
240 MS=1024 

250 REDIM X(MS) ,Y(MS,2) 

260 P4i$="VOLTAGE” 

270 P2$="DIGDE NUMBER" 

280 P3$="VOLTAGE VS. DIODE" 
290 P4$="PLTDTA RESULTS" 
300 ! 

310 ASSIGN#® 4 TO Q3$ 

320 FOR I=4 TO MS-7 STEP 9 
330 READ# 4 ; A,B,C,D,E, 
340 Y(I,4)=A @ YCI+4,4)= 
(I+4,14)=E£ @ Y(I+5,4)=F 
350 Y(1+6,4)=G @ YCI+7, 4) =H 

360 Y(I,2)=R @ Y(I+i,2)=S @ Y(1+2,2)=T @ Y(1+3,2)=U @ Y 
(I1+4,2)=Y @ Y(I+S,2)=W 

370 Y(I+6,2)=X @ Y(1+7,2)=Z 

380 NEXT I 

390 M7=MAXAB(Y) 

400 DISP "DATA INPUT COMPLETE" 

440 LIMIT 0,250,0,1980 

420 FRAME 

430 LOCATE 20,120,20,80 

440 FRAME 


F,G,H 
B @ Y 3,4)=D @ Y 


100 





450 SCALE 0,MS,-.1,m7+.4 

460 CSIZE 3,.5,0 

470 DEG 

480 LDIR @ 

490 FXD 0,2 

S00 LAXES MS/20,.02,0,0,2,5,5 
S40 MOVE .6SXMS, .93K(M7+. 4) 

S20 LABEL P4$ 

530 MOVE .6SKMS,.84k(M7+.4) 

S40 LABEL P3$ 

SSO LABEL Tis 

S60 MOVE MS/2,-. 055 

S70 LABEL P2$ 

S80 LDIR 90 

S90 MOVE ~(.4%MS), .SXM7 

600 LABEL Pis 

610 LDIR 0 

620 M=4 

630 GOSUB 760 ! PLOT FIRST ARRAY 
640 DISP " FIRST ARRAY PLOTTED!" 
650 M=2 

660 GOSUB 7460 ! PLOT SECOND ARRA492 MOVE .465%MS, .74&(M7 
+. 4) 

670 DISP “ SECOND ARRAY PLOTTED!" 
680 DISP “ PLOT COMPLETE! “ 

690 DISP "DO YOU WISH TO PLOT THE OTHER PHOTO ARRAY FIL 
E? (Y/N) " 

700 INPUT N$ 

740 IF N$="Y" THEN GOTO 410 

720 CLEAR 

730 DISP "LOAD AND RUN (REDUCE) TO PROCEED WITH DATA RE 
DUCTIONS., * 

750 END 

760 ! PLOT DATA 

770 X(4)=4 

780 MOVE X(4),Y¢4,4) 

790 FOR I=i TO MS STEP J 

800 X(I)=I 

810 PLOT X(1I),Y(1,M) 

820 NEXT I 

830 PENUP 

840 RETURN 
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po! REDUCE 


20 CLEAR 7 
30 DISP “REDUCE LOADED" 
40 ! 


SO ! BERT HANSEN 

bO ! THIS PROGRAM REDUCES THE ARRAY DATA 

70 OPTION BASE i 

80 DIM A(2,2),B(2),Y4(2) 

90 COM Q3$(121,Q2$(121,T1$(201],Di$li2],D2$(4121,4 

100 QUTPUT 709 ,"TD" 

440 ENTER 709 ; TAs 

120 CLEAR 

i30 DISP "DATA REDUCTION BEGINS" 

140 DISP “ ENTER THE INDEX OF REFRACTION=?(AIR M=1.0, W 

ATER M=i.35)" 

4S0 INPUT M 

160 DISP “INSTALL PAPER AND PEN ON PLOTTER TO PLOT I() 
vs " 

{70 DISP “PRESS[CONT] WHEN DONE" 

180 PAUSE 

190 DISP "ENTER THE FILE NAME OF DATA TO BE REDUCED." 
200 DISP “INCLUDE THE DISC LOCATION ON FILENAME ( NAME: 

D701)" 

240 INPUT A3 

220 D9$=":D704" 

230 L=.000000488 ! WAVELENGTH IN METERS 

240 SHORT T4(1024) ,14(4024) ,12(4024) 

250 MS=1024 ! CURRENT NUMBER OF ACTUAL READINGS 

260 D8=1.5365 |! MMETERS 

270 F8=S26 |! MMETERS 
280 Di=.025 ! MMETERS 

285 GCLEAR 

290 PLOTTER IS 70S 

300 LIMIT 0,250,0,180 

340 LOCATE 20,120,20,80 

320 FRAME 

330 SCALE 0,.04,-.4,4 
340 CSIZE 3,.5,0 

350 DEG 

360 LDIR 4 

370 FXD 3,4 

380 LAXES .002,.4,0,0,2,2,5 

390 DISP "DO YOU WANT TO PLOT A COMPARISON CURVE ( UNIC 
UR)? (Y/N)" 

400 INPUT NS 

440 IF N$="Y" THEN GOSUB 4440 

420 ASSIGN#® 4 TO AS 

430 FOR I=4 TO MS-7 STEP 8 

440 READ 4; A 

450 [2(1)=A @ I 





4)=E @ L2Q(I+S)=F @ I2(1+6)=G6 @ [2¢1+7)=H 
460 LAC L)=R @ LTACI+i)=S @ LiCI+2)=T @ LA(I+3)=U @ Licr+ 
4)=V @ Ti(I+S)=W @ Li CI +6)=X @ 11 ¢1+7)=Z 
470 NEXT I 

480 MAT I4=114-12 

490 ZI=AMAX(14) 

S00 Z0=AMAXROW 

S40 PRINT "MAXIMUM OUTPUT AT DIODE#=";20 
S20 DISP "ENTER THE INCREMENT AFTER PEAK (J=#START,K= 
STOP FIT)” 

S30 INPUT J,K 

S40 NO=K-J 

5S0 Xi=0 @ X2=0 @ X3=0 @ Yi=0 @ Y2=0 

S60 FOR N=J TO K 

S70 Ti (N)=N+61.46 

S80 X1=Xi+Ti(N) 

S90 X2=X2+TL(N)*2 

600 XB=X3+Ti (ND KILN) 

640 YI=YAt+1Ii(N) 

620 Y2=Y2+Ii(N)*2 

630 NEXT N 

640 ACL, 1)=ND 

650 A(41,2)=X4 

660 A(2,1)=X4 

670 A(2,2)=X2 

680 Yi(4)=Y4 

690 Y4(2)=X3 

700 MAT B=SYS(A,Y4) 

740 PRINT “DELTA="3B(4) 

720 PRINT “SLOPE="35(2) 

730 AB=B(4) 

740 R2=(NOKX3B-XiKY4)*2/¢ (NOKX2-X4°2) KC NOKY2-Y442)) |! CO 
RRELATION 

750 PRINT 

760 PRINT "CORRELATION R*2=";3R2 

770 MAT I4=(4/A3) X14 

780 ! I4 ARRAY IS THE FINAL PLOTTING NORMALIZED ARRAY 
790 ! NOW PLOT Ii VS, Ti 

800 FOR I=i TO MS 

840 Ti(I)=(1*,025+D8)/(F8*M) |! THETA 

820 PLOT T4i(I),14¢1) 

830 NEXT I 

840 PENUP 

850 MOVE 0,0 

860 MOVE .02,-.3 

870 LABEL "THETA (RADIANS) " 

880 MOVE -.004,.2 ® DEG @ LDIR 90 

890 LABEL "NORMALIZED INTENSITY" 

900 LDIR 0 @ MOVE .015,.8 

940 LABEL “INTENSITY VS. THETA" 
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920 
930 
940 
950 
760 
970 
9980 
990 


MOVE 


LAKEL 


MOVE 


LABEL 


,045,.9 


eRe DUCES RESULTS" 


,045,.7 


Tis 


QS$="REDUCS" 


DISP 


IS THIS A CALIBRATION RUN? (Y/N)" 


INPUT N$ 
IF NS="N" THEN GOTO 1010 


1000 
1010 


QS$="CALTE" 
DISP "ENTER THE NAME OF THE DATA FILE TO STORE THE 


DATA CINTENSITY VS,.THETA)" 


1020 DISP “INCLUDE THE DISC LOCATION ON FILENAME ( NAME 
:D704)" 

1030 INPUT Q33 

1040 CREATE Q3$,1024, 16 

1050 PRINT “INTENSITY VS.THETA STORED IN FILE=";Q3% 
{060 ASSIGN# 2 TO Q3$ 

1070 FOR I=4 TO MS 

1080 PRINT# 2 ; Ti(I),1icd) 

{090 NEXT I 

41400 ASSIGN# 2 TO 3 

1440 DISP "DATA STORED ON DISK IN FILE -- ";Q3$ 
1420 CHAIN QS$ 

1430 END 

4440 ! COMPARE WITH UNICUR 

11450 ASSIGN# 2 TO "“UNICUR" 

{160 FOR I=it TO 4100 STEP 2 

1170 READ# 2 3 Ti(I),I4¢1),TACI+4) , Li C144) 

1480 NEXT I 

1490 DISP "ENTER EXPECTED VALUE OF Dbar (23.98E-6)ETC.”* 
4200 INPUT D 

1210 FOR I=i TO 99 


1220 
1230 
1240 
1250 
1260 


TACT) =TACIDKL/CPIRDD 
PLOT Ti(I),14(1) 
NEXT I 

PENUP 

RETURN 
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10 

20 

30 

40 

50 

&0 

70 

80 

90 

100 
110 
120 
130 
140 
5) * 
150 
160 
05) 
170 
180 
190 
200 
210 
220 
2.30 
240 
250 
260 
270 
280 
290 
500 
310 
320 
330 
3490 


1 CALIR 
OPTION KASE 1 
COM Q3$(4121,Q2$(121,71$(20),D4$l12) ,D2$12],m 


D9$=":D704" 
1 CALIBRATION PROCEDURE 

1 BERT HANSEN 
CLEAR 
GCLEAR 
DISP "ARE YOU RUNNING CALIB SEPARATELY? (Y/N)" 

INPUT N$ 

IF N$="N" THEN 460 

DISP “ENTER DATAFILE NAME AND DRIVE ADDRESS" 

INPUT Q3$ 

DISP “ENTER INDEX OF REFRACTION (AIR=1.0 ,WATER=14.3 
INPUT M 

DISP "ENTER THE ADDRESS OF PLOTTER. (CRT=1, 722SB=7 
INPUT P 

PLOTTER IS P 

IF P=70S THEN LIMIT 0,250,0,180 

LOCATE 20,120,20,80 

FRAME 

SCALE 0,10,-.1,4 

LAXES .S,.2,0,0,2,5,% 

} 

DIM ACi4,414),BC44) ,WC20) ,YiCit) ,Y2C44) 

SHORT X(41024) ,Y(1024) 

MAT A=ZER 

MAT B=ZER 

MAT Yi=ZER 

MAT W=ZER 

MS=41024 

DISP " INDEX OF REFRACTION IS ="5M 

DISP "USING ¢ ";Q3%;" >" 

DISP "A POLYNOMIAL REGRESSION ROUTINE OF (J)TH ORDE 


R WILL BE PERFORMED. " 


4516 


DISP "ENTER ORDER OF POLYNOMIAL( REASONABLE), " 


360 INPUT Ji 


370 


ASSIGN# 1 TO Q3% 


380 FOR I=4 TO MS STEP 2 
390 READ# 1 5 X(I),YCI),X(L+1), YC 144) 


400 NEXT I | 

410 DISP "ENTER Dbar TO USE IN DATA REDUCTION (e.g. 23 
_7@e-6)" 

420 INPUT D3 


430 L=.000000488 


440 


F=526 ! FOCAL LENGTH 


450 K6=D3KePT/L 


460 


FOR T=i TO MS 
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470 
480 
490 
$00 
510 
520 
5.30 


550 
560 


pas, 


570 
$80 
$90 
600 
610 
620 
630 
640 
650 
660 
6790 
680 
690 
700 
710 
720 
730 
740 
7539 
7690 
7790 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
710 
920 
930 
940 
950 
960 


X(1)=X (1) *K6 

Peat yy «1 ) 

NEXT I 

' I YS.THETA BAR 

PENUP 

DISP "ENTER THE MINIMUM THETA BAR OF INTEREST." 
INPUT J 

J=MKI/ (236KD3)-614.46 

K=4,5365 

DISP “ENTER THE MAXIMUM THETA BAR OF INTEREST .(ba 
Ss)" 

INPUT M4 

Mi=(MI¥LKFKM/ (PIXD3)-K)/, 02533 

PRINT “MIN DIODE ="3;J 

PRINT “MAX DIODE ="3M4 

J2=Ji+i 

FOR I=4 TO 24J4 

j 


S0=0 

FOR M=J TO Mi 

f 

SO=X(M)*I+S0 

NEXT M 

WCT)=S0 

NEXT I 

1 

FOR I=4 TO J2 

S0=0 

FOR N=J TO Mt 

SO=SO+Y (N)KX(N) *( I-14) 
NEXT N 

Yi¢1)=S0 

NEXT I 

NO=Mi-J ! # OF DATA PAIRS 
{ 

REDIM ACJ2,J2),B(J2) ,WC2xT4) ,Y4(I2), Y2(T2) 
A(4,4)=NO 

FOR I=2 TO J2 ! FIRST ROW 
ACL, 1)=W(I-4) 

NEXT I 

| 


FOR I=2 TO J2 ! FIRST COLUMN 
ACL, 4)=WC(I-4) 

NEXT I 

j 


M=0 

FOR I=2 TO J2 ! CORE 10*19 
FOR K=2 TO Je 
ACT, KO =WKK +) 

NEXT K 

M=M+1 
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v7 OO NEXT 1 , 
780 GOSUB 1300 ! 


990 

1000 
1040 
4020 
1030 
1040 
4050 
1060 
10790 
1080 
10990 
1100 
14140 
1120 
1130 
1140 
1150 
1160 
1170 
£180 
1490 
4200 
4210 
1220 
4230 
1240 
4250 
1260 
1270 
1280 
2290 
4300 
43140 
1320 
1330 
1340 
1356 
4360 
4370 
1380 


CORRELATION 


MAT K=SYS(A,Y4) 
PRINT 


PRINT 
PRINT 


"POLYNOMIAL COEFFICIENTS=; " 


MAT PRINT 


° 


! COMPUTE THE CALIBRATION CURVE 


O=.4 


FOR I=i1 TO 100 
XCI)=1%$ 


50=0 


FOR N=14 TO Je 
S0=SO+BCN) KXCT)*(N-14) 
NEXT N 


Y(I)= 
PLOT 


S0 
Lh yo 


NEXT I 


DISP 
DISP 


PVEGMPLETE !" 
"ENTER THE DESIRED FILE NAME TO STORE DATA." 


INPUT Q3$ 

CREATE Q3$,100,46 

ASSIGN# 2 TO Q3$ 

! STORE SMOOTH CAL CURVE ON DISC. 
FOR I=4 TO 100 STEP 2 


PRINT# 2 ; 


AOL) PGineeA 6 Peal +1 ) 


NEXT I 


DISP 


"DATA STORED ON DISC" 


ASSIGN# e@ TO xX 


CHAIN 


END 


"REDUC2" 


! COMPUTE CORRELATION FACTOR 


DISP 


“ CORRELATION " 


SHORT X4(4024) ,Y4(4024) 
MAT X4=ZER 

MAT Y4=ZER 

MAT X4=X( 3:64) 

MAT Y4=Y(I:M4) 
C1=DOT(X4,X4) 
C2=DOT(Y4,Y4) 


LO 








1390 C3=DOT(X4,Y4) 
1400 C4=Ci*2 

1440 CS=C2*2 

1420 C6=SUM(X4) 

{430 C7=SUM(Y4) 

1440 N2=N0*C3-C&x*C7 

{450 D2=SQR¢{ (NOKC1-CH*2) K(NOKC2-C7*2)>d 
1460 R=N2/D2 ! CORRELATION FACTOR 

{470 R2=R*2 

4480 PRINT " CORRELATION FACTOR” 

1490 PRINT “R‘42=";R2 

1500 RETURN 
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10 !§ REDUC2 

20 ! 

30 ! BERT HANSEN 

40 ! THIS PROGRAM WILL RETRETIVE THE UNIVERSAL CURVE DA 


TA FROM THE DISC. 

SO ! IT WILL THEN PLOT THE CURVE AND THEN RETRELVE AND 
PLOT THE CALIBRATION CURVES. 
60 | 

70 CLEAR 

80 OPTION BASE 4 

90 COM Q3$(412],Q2$(421,T1$(20) ,Di$l12],D2$1121,m 
400 D9$=";D704" 

1140 L=.000000488 

120 F=S526 

130 Ki=1.5365 

140 PLOTTER LS 705 

iS0 CSIZE 3,.4,0 

160 QUTPUT 709 ;"TD" 

170 ENTER 709 ; Tis 

480 SHORT X(100),Y(4100) 

190 SHORT T(1024) ,1(41024) 

200 P4$="THETA BAR“ 

240 PS$="INTENSITY" 

220 DIM P3$[26] 

230 P3$="INTENSITY VS. THETA BAR" 
240 LIMIT 0,250,0,180 

260 LOCATE 20,4120,20,80 

270 FRAME 

280 SCALE 0,10,0,3 

290 LINETYPE i 

300 E=4 

340 E3=3 

320 E1$="A,.3D" @ E2$=".3D,A" 
330 Y=4 @ Y2=10 

340 GOSUB 1230 ! LABEL 

350 GOSUB 1300 ! LOG SCALE 
360 E1$="A,.2D" @ E2$=".2D,A" 
370 E=10%*E 

380 Y=20 @ Y2=100 

390 COSUB 1300 

400 Ei$="A,D.D" @ E2$="D.D,A" 
440 Y=200 @ Y2=1000 

420 E=10*E 

430 GOSUE 12308 

440 DEG 

450 LDIR -99 

460 FOR X=0 TO 10 

470 ! YAXIS x 

480 MOVE X,-.4 

490 LABEL USING "DD" ; x 


LOS 





500 
510 
oe 0 
530 
540 
So) 
560 
579 
580 
590 
600 
610 
62 

& 30 
640 
650 
660 
§70 
680 
699 
7090 
710 
720 
730 
740 
750 
769 
770 
780 
770 
800 
810 


NEXT xX 

LDIR Q 

DISP "SCALE COMPLETE" 
SCALE 0,10,6,14 

CSIZE 4,.5,0 

FXD 4,2 

MOVE S,-.45 

LABEL P4$ 

DEG 

LDIR 90 

MOVE -1,.5 

LABEL PS$ 

LDIR 0 

MOVE 7,.85 

| LABEL P3$ 

MOVE 7,.9S 

' LABEL "REDUC2 RESULTS" 
MOVE 7,.75 

|! LABEL Ti$ 
QS$="UNICUR" 

DISP " LABELING COMPLETE" 
DISP "CHANGE PENS AND [CONT)" 
PAUSE 

SCALE 0,10,-3,9 
MS=1024 

ASSIGN#® 4 TO Qss 

FOR I=i TO 99 

READS 4 5 X(1),YC(I) 
Peon XCL) ,LGT(Y<1)) 
NEXT I 

PENUP 

ASSIGN® 4 TO x 


820 DISP " DO YOU WANT TO PLOT A CALIBRATION CURVE? (Y/ 
N)* 

830 INPUT NS 

340 IF N$="N" THEN GOSUB 930 

850 DISP " ENTER THE FILE NAME OF THE CALIBRATION CURVE 
TO BE PLOTTED. (CALi,ETC)" 

B60 DISP “ INCLUDE THE DISC DRIVE WHERE CAL FILE IS STO 
RED--E.G. -NAME: D704" 

870 INPUT Q3% 

B80 ! DISP "ENTER THE NUMBER OF DATA POINTS IN THE FILE 
.(400,4024,ETC)" 

890 ! INPUT M2 ! CAL FILES HAVE 100 ELEMENTS 

900 M2=400 

940 GOSUB 980 ! CALC THETA BAR AND PLOT |! 

920 GOTO 820 

930 LINETYPE 4 

940 CLEAR 
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750 
760 
970 
980 
290 


1000 
10410 
1020 
10350 


TEREST. 


1040 
1050 
1060 
1070 
1080 
1090 
1100 
111090 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
4190 
1200 
{210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
4320 
1330 
1340 
1350 
1360 


DS 
DISP 
END 


DISP 
L RUN. 


Pree to COMPEER TE. b" 
"LOAD REDUCS TO COMPUTE D32 FROM A DATA SET." 


READ AND PLOT CAL 

“ENTER THE INDEX OF REFRACTION USED FOR THE CA 
(AIR=4.0, WATER=1.35)" 

INPUT M 

DISP “ENTER THE Dbar USED TO COMPUTE THETA KAR.” 
INPUT D 

DISP “ENTER THE DIODE NUMBERS FOR THE REGION OF IN 
(P4,P2) TYPICALLY (50,400) " 

INPUT P1,P2 

H9=DxP1*,025/(L¥*F MD 
P4=H9K(64.46+P 4) 

P2=H9K(64,.46+P2) 

PRINT "LOWER LIMIT OF THETA BAR=";P4 
PRINT 

PRINT "UPPER LIMIT OF THETA BAR=";P2 
ASSIGN#® 4 TO Q3$ 

FOR I=4 TO M2 

READ@ 4 ; TC(I),ICDD 

NEXT I 

MOVE T(4),1¢4) 

FOR I=4 TO 100 

IF T(I)<P4 THEN 4200 

IF T(I)>P2 THEN 1200 

Meet) 90) THEN PLOT 1¢(1),LGT(1IC1)> 
NEXT I 

PENUP 

RETURN 

! LABEL 

Y3=Y/41000 

MOVE 10.4,LGT(Y) @ LORG 2 

LABEL USING Ei$ ; "-",Y3 

MOVE ~.4,LGT(Y) @ LORG @ 

LABEL USING E2$ ; Y3,"-' 

RETURN 

FOR Y=Y TO Y2 STEP E 

E2=LGT(Y) 

1 XAXIS LGTCY) 

IF RMDCY,40%E)#0 THEN 4350 

GOSUB 1230 ! LABEL 

NEXT Y 

RETURN 
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10 ! REDUCZ 

40 | 

SO ! THIS ROUTINE COMPUTE THE VALUE OF D32 FOR THE DA 
TA IN THE FILE "Q$" . 

60 ! BERT HANSEN 43 JULY 82 

70 | 

80 OPTION BASE i 

90 MS=1024 

400 COM Q3$l12] ,Q2$(421,74$(201] ,D4i$(421,D2$([14121,m 

140 CLEAR 

120 D9$=":D704" 

430 SHORT R(20),X4¢100),Y4(100) 

440 SHORT X(4024) ,Y(41024) ,D3(20) ,T3(20) ,T(4024) 

150 DISP “ENTER THE FILENAME OF THE DATA TO BE REDUCED 

(DATAL:D704,ETC)" 

460 INPUT Q3$ 

170 ASSIGN# 2 TO Q34 

180 FOR I=1 TO MS-3 STEP 4 

190 READ@® 2 ; X(I),YCL) ,XCL+4),YCI44) ,XCI+2) ,YCI42) ,XC1 
#3) ,Y(1+3) 

200 NEXT I 

240 ! FOR T=4 TO 1024 

220 1 X(I)=X(1)/1998. Sa 

230 ! NEXT I 

240 ! THETA AND INTENSITY ENTERRED 

250 L=,000000488 

260 Ki=L/PI 

270 DISP “ENTER THE NAME OF THE DATA FILE OF THE REFERE 
NCE ARRAY. (UNICUR ,CAL4,CAL2,CAL3)" 

280 INPUT Qis$ 

290 M4=4100 

300 ASSIGN# 4 TO ais 

340 FOR I=4 TO M4-4 STEP 2 

320 READS 4 5 X4(1),YiCIT) ,XiCI +4) ,YACI +4) 

330 NEXT I 

340 DATA .5,.3,.2,.4,.09,.08,.07,.06,.05,.04,.03,.02,.0 
4,.008,.00S 

350 FOR N=1 TO iS 

360 READ R(N) ! REF INTENSITIES 

370 NEXT N 

380 DISP “ENTER THE NUMBER OF DIODES TO BE OMMITTED CLO 
SE TO THE CENTERLINE. " 

390 INPUT J 

400 ! D32 CALCULATION 

440 Ji=2 

420 FOR D=4 TO 4S 

430 Rt=R(D) ! REF INTENSITY 

440 OL1=AMIN(Y) |! SMALLEST INTENSITY 

450 IF R4<04 THEN GOTO 1250 

460 FOR I=Ji TO 100 


2 





470 
480 
490 
560 
510 
520 
530 
540 
550 
560 
570 
$80 
590 
600 
610 
620 


IF Yi(I))R4i THEN GOTO 610 

IF YiCL1))=Yi(1-4£) THEN GOTO 600 
| FIND THETA EAR 
H2=Yi¢1I-4)-YiCL) 

K2=X4(1)-Xi (I-14) 

H3=Y4i(I-41)-R4 

K3=K 2% (H3/H2) 

T3(D)=X4¢1-4)+K3 ! THETA BAR 
PRINT “I="sI;YiCI) jR45D 

PRINT "L="sI;T3(D) ;X4¢1-4)4+K3 
PRINT 

Ji=I 

GOTO 620 

PRINT " ZERO OR POS SLOPE" 

NEXT I 

! FIND THE INTENSITY CORRESSPONDING TO DESIRED INT 


ENSITY 


630 
640 
650 
569 
670 
680 
690 
700 
7190 
720 
730 
740 
7590 
760 


779 
780 
790 
800 
819 
820 
839 
840 
850 
B60 
870 
880 
8990 
900 
710 
920 
930 
940 
950 


M9=1000 

L9=4 

FOR I=J TQ M9 

IF Y¥(I)>0 THEN 696 

PRINT “INTENSITY VECTOR GOING NEGATIVE AT I=";1 
PRINT “INCREASING J=J+5" @ J=J+5 ® GOTO 639 
Z=0 

FOR M=I-S TO [+5 

L=Z+Y(M) 

NEXT 

Z=,1X2 

TF Z>Ri THEN GOTO 790 

! 21S LESS THAN Ri 

J=I ! STARTING POS FOR NEXT SEARCH 


DISP “ I FOUND IT" 

GOTO 810 

L9=Z ! LAST AVERAGE: 

NEXT I 

1 FIND THETA 

K7=X(J)-X( I-14) 

H7=L9-Z 

H6=L9-Ri 

K6=K7*(H6/H7) 

T4=X(J-4)+K6 ! THETA 

! COMPUTE D32 
D3(D)=KiXT3(D)/T4 ! D32 

PRINT 

PRINT "FOR INTENSITY=";R(D) 
PRINT "THETA= ";T4;" RADIANS" 
PRINT "THETA BAR= "3; T3C(D) 
PRINT "D32= "3;D3(D)*1000000;" MICRONS” 
NEXT D 

DS=0 


ies 





960 FOR Ti=i TO D 
970 DS=DS+D3(T1) 
980 NEXT [14 

990 DS=DS/Dx*i 000000 


{000 PRINT “D32 AVERAGE=";D5;" MICRONS “ 

4010 PRINT 

{020 DISP "DO YOU WANT TO COMPARE DATA TO ANOTHER REFER 
ENCE FILE? (Y/N)" 

1030 INPUT N$ 

4040 IF NS="N" THEN 1080 

1050 DISP “ENTER THE FILENAME OF THE REF. DATAFILE (CAL 
1A,ETC)" 

1060 INPUT Q1$ 

1070 GOTO 290 

1080 DISP "DO YOU WANT TO REDUCE ANOTHER DATAFILE? (Y/N 
a 

4090 INPUT NS 

{400 IF NS="N" THEN 4480 

1440 DISP “ENTER THE FILENAME OF THE DATAFILE TO BE RE 
DUCED. (DATAL,ETC)" 

4420 INPUT Qs 

1430 ASSIGN# 1 TO Q3 

1440 FOR I=4 TO 41024 STEP 4 

£1450 READ#® 4 3 X(I),YCI),XCT+4),YC1#4) ,XC14+2) , YC I+2) , X¢ 
+3) ,Y(I+3) 

{460 NEXT I 

1470 GOTO 400 ! REPEAT D32 CALCULATIONS 

1480 DISP “REDUCTIONS COMPLETED" 

1490 DISP “WANT TO DO ANOTHER RUN?” 

{200 INPUT N$ 

1240 IF N$="N" THEN 4230 

1220 GOTO 10 

{230 CHAIN “REDUC4" 

4240 END 

1250 ! REF INTENSITY IS LESS THEN SMALLEST INTENSITY 
{260 ! DON’T SEARCH JUST COMPUTE D32 AVERAGE TO THIS PO 
INT 

{270 DS=0 

1280 FOR Ii=i TO D-i 

1290 DS=DS+D3( 11) 

{300 NEXT If 

1340 DS=DS/(D-4) ! D32 AVERAGE 

{320 GOTO 1000 
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i 0 
20 
30 
40 
TA 
50 


60 


REDUC4 


BERT HANSEN 
THIS PROGRAM WILL RETREIVE THE UNIVERSAL CURVE DA 


! 
! 
| 
FROM THE DISC. 
i 


TY WILL THEN PLOT THE UNIVERSAL CURVE. 


! THE USER MAY THEN PLOT INTENSITY VERSUS THETA HAR 


WITH D32 ENTERRED FROM KEYBOARD. 


79 


° 


80 CLEAR 
90 GCLEAR 


100 
i190 
120 
130 
140 
150 
160 


{70 
180 
190 
200 
2id 


aed 


230 
240 
250 
260 
270 
280 
290 
300 
340 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 


460) 


OPTION BASE if 
COM Q3$(412),Q2S012)],7T1$020),D1$li2] , D2stie2l,M 
D9$=";D704" 
=,000000488 
F=526 
Ki=41.5365 
DISP “ENTER THE PLOTTER ADDRESS (CRT=1, 722S5B=705)" 


INPUT P 

PLOTTER IS P 

IF P=70S THEN LIMIT 0,250,0,180 
CSIZE 3,.4,0 

OUTPUT 709 ;"TD" 

ENTER 709 3 Tig 

SHORT X(400),¥¢(400) 

SHORT 17(1024) ,1(4024) 
P4$="THETA BAR" 
PS$="INTENSITY" 

DIM P3$(26] 
P3$="INTENSITY VS. THETA BAR” 
FRAME 

LOCATE 20,120,20,380 

FRAME 

SCALE 0,10,0,3 

LINETYPE 14 

E=4 

E3=] 

Ei$="A,.3D" @ £2$=",.3D,A" 
Y=4 @ Y2=40 

GOSUB 1260 ! LABEL 

GOSUB 4330 ! LOG SCALE 
Ei$="A,.2D" @ E2$=",2D,A" 
E=410xE 

Y=20 @ Y2=400 

GOSUER 4330 

ELisg="A,D.D" @ EF. ="D.D,A" 
Y=200 @ Y2=1000 

E=40*E 


a1 





470 GOSUB 1330 

480 DEG 

490 LDIR -99 

S00 FOR X=0 TO 10 

S40 YAXIS X 

S20 MOVE X,-.4 

S30 LABEL USING "DD" ; X 

S40 NEXT X 

SSO LDIR ¢ 

S60 DISP "SCALE COMPLETE“ 

S70 SCALE 0,10,0,14 

S80 CSIZE 4,.S5,0 

S90 FXD 4,2 

600 MOVE S,-.45 

640 LABEL P4$ 

620 DEG 

630 LDIR 99 

640 MOVE -1,.S 

650 LABEL PS$ 

660 LDIR 0 

670 MOVE 7,.985 

680 LABEL P3% 

690 MOVE 7,.95 

700 LABEL "REDUC4 RESULTS" 

740 MOVE 7,.75 

720 LABEL Tis 

730 QS5$="UNICUR" 

740 DISP “ LABELING COMPLETE" 

750 DISP “CHANGE PENS AND [CONT]" 
760 PAUSE 

770 SCALE 0,10,-3,0 

780 MS=1024 

790 ASSIGN# { TO QS$ 

800 FOR I=4 T0 99 

840 READ@ 4 ; X(1),YCI) 

820 PLOT X(I),LGT(Y(I)) 

830 NEXT I 

840 PENUP 

850 ASSIGN# 1 TO x 

860 DISP " DO YOU WANT TO PLOT A DATA FILE CURVE? (Y/ 
N>“ 

870 INPUT N$ 

880 IF N$="N" THEN GOSUB 970 

890 DISP " ENTER THE FILE NAME OF THE DATA FILE TO HE 
PLOTTED. (NAME;D704)" 

900 DISP " INCLUDE THE DISC DRIVE WHERE CAL FILE IS STO 
RED--E.G. -NAME:D7014" 

940 INPUT Q3¢ 

920 M2=1024 

930 GOSUB 10410 ! CALC THETA BAR AND PLOT ! 


SG 





940 

9590 

260 

TAY 

980 

790 

4000 
10410 
1020 
AL R 
1030 
1040 
BAR . 
10590 
4060 
TERE 
10790 
1080 
10990 
1400 
1410 
1120 
14130 
1140 
1150 
1160 
4170 
T+3) 
1180 
1199 
1200 
4210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
£290 
4300 
13190 
1320 
1330 
4340 
1350 
1360 
4370 
1380 
13990 


IF P=705 THEN 9?70 

GRAPH 

COPY 

BENETYPE 1 

CLEAR 

DISP 
END 
|! READ AND PLOT CAL 


feo COMPLETE | 


DISP "ENTER THE INDEX OF REFRACTION USED FOR THE C 
UN. (AIR=41.0, WATER=1.35)" 

INPUT ™ 

DISP “ENTER THE D32 TO BE USED IN COMPUTING THETA 
INPUT D 

DISP “ENTER THE DIODE NUMBERS FOR THE REGION OF IN 
ST. (Pi,P2) TYPICALLY (50,400)" 


INPUT P1i,P2 

H9=DxP Ik, 025/ (LF KM) 

P3Z=H9*(641.46+P4) 

P4=H9X(64,45+P2) 

K6=D¥kPI/L 

PRINT “LOWER LIMIT OF THETA BAR=";P3 

PRINT 

PRINT “UPPER LIMIT OF THETA BAR=";P4 

ASSIGN® 4 TO Q3% 

FOR I=i4 TO M2 STEP 4 

READS 4 » TCL), ICL), TCI4+£i) ,LCI+i) ,TCI+2) ,L¢1 +2) , TC 
»1¢1+3) 

NEXT I 

MOVE T(4),I(4) 

FOR I=Pi TO P2 

TCL) =K64TCI) 

IF I1¢(1))0 THEN PLOT T(I),LGT(I(I)) 

NEXT I 

PENUP 

RETURN 

1 LABEL 

Y3=Y/41000 

MOVE 10.4,LGT(Y) @ LORG 
LABEL USING Ei$ 3; "-",Y3 
MOVE -.4,LGT(Y) @ LORG 8 
LABEL USING E2$ ; Y3,"-" 
RETURN 

FOR Y=Y TO 
E2=LGT(Y) 
XAXIS LGTCY) 

IF RMD(Y,10*E)#0 THEN 1380 
GOSUB 4260 ! LABEL 

NEXT Y 

RETURN 


> 


2 


® 
? 


Y2 STEP E 
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